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Abstract. Let T> denote the class of bounded real analytic plane domains 
with the symmetry of an ellipse. Under generic conditions, we prove that if 
^li^2 £ "D and if the Dirichlet spectra coincide, Spec(Oi) = Spec(p,2), then 
f2i = Q.2 up to rigid motion. 
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1. Introduction 

In this paper we give a positive solution to the inverse spectral problem for the 
class T> of analytic axi-symmetric plane domains satisfying: 

• S7 is real analytic 

• f2 is Z2 x Z2 — symmetric 

• at least one axis of is a non — degenerate bouncing ball orbit 7 

• L~ i has multiplicity one in Lsp(Jl) . 

Let Spec(Jl) denote the spectrum of the Laplacian Aq of the domain with Dirich- 
let boundary conditions. Our main result is the following theorem (announced in 
|H): 

Theorem 1.1. Spec: V i-> is 1-1. 

Let us clarify the assumptions. The symmetry assumption is that there is an 
up/down reflection symmetry across a horizontal axis and a left/right reflection 
symmetry across a vertical axis. Both axes intersect the boundary at right angles, 
hence are projections to SI of 'bouncing ball orbits' of the the billiard flow G l on 
T*fi, with the usual law of reflection at the boundary. Associated to any periodic 
reflecting ray f3 of G* is its Poincare map Vp, defined as the first return map on 
a local transversal to f3 in S*fl, and its linear Poincare map Pp — dp/QsVp. A 
periodic orbit (3 is said to be non-degenerate if no eigenvalue of Pp is a root of 
unity. Thus, it is non-degenerate elliptic if the eigenvalues of Pp have the form 
{e ±ia } with a/ir ^ Q, or non-degenerate hyperbolic if they have the form {e ±A } 
with A 7^ 0. Our non-degeneracy assumption is that at least one of the axes, which 
we will denote by AB, is the projection to Jl of a non-degenerate bouncing ball orbit 
7. Also, Lsp(Sl) denotes the length spectrum of f2, i.e. the set of lengths of periodic 
billiard trajectories (including the boundary). Our length spectrum assumption is 
that 7 is of multiplicity one in Lsp(£l) 7 i.e. that it is the unique trajectory of its 
length L 1 . With no loss of generality we will assume AB is the vertical axis, and 
will denote its length by L; thus L 1 = 2L. 
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The proof of Theorem (1.1) is based on the method of normal forms, which was 



introduced into inverse spectral theory by Colin de Verdiere [CV| and by Guillemin 
|G|, The basic idea of Guillemin jG| is to construct a quantum analogue of the 
Birkhoff normal form of Aq around each closed geodesic and to prove that the 
coefficients B^ of the normal form are spectral invariants. The latter is accom- 
plished by relating these coefficients B^ to the wave invariants o 7 fc of An at 7, 
i.e. the coefficients of singularity expansion of the trace of the wave group U(t) 
at the time t = L~ i . In Q, Guillemin constructed the normal form and proved 
that it is a spectral invariant in the case of non-degenerate elliptic closed geodesies 
on boundaryless manifolds. A somewhat different construction of the normal form 



and a new proof that the coefficients are spectral invariants was given in |Z.l] [ p5.2| 
for general non-degenerate closed geodesies. The inverse spectral problem is then 
reduced to determining the metric from the normal form. 

The latter inverse problem remains difficult since the different closed geodesies 
do not easily 'communicate' with each other and since the B 1 k , s for a fixed 7 (and 
its iterates) do not appear to give enough information to determine the metric, 
even locally. Therefore it is natural to consider the problem first for classes of real 
analytic metrics with one functional degree of freedom, where it is plausible that the 
normal form coefficients at just one closed geodesic should determine the metric. 
This motivated our study in |p5.3|| of real analytic 'simple' surfaces of revolution. 
We proved there that a simple analytic surface of revolution is determined by the 
normal form of its An along the 'meridian torus' of closed geodesies. For further 



background on the method of normal forms we refer to the expository artice [Z.4|. 

In this paper, we extend the method of normal forms to the case of bounded 
analytic plane domains with a non-degenerate bouncing ball orbit 7. We do not 
use the symmetry assumption on the domain in the construction of the normal form 
at 7, but only in the last step of deducing the domain from the normal form. Since 
the construction of the normal form and the various technical pitfalls may make 
the proof difficult to follow, let us give a brief summary here of the main ideas. See 



also [Z.5] for a somewhat less technical outline. 

As mentioned above, the main idea is to introduce a notion of Birkhoff normal 
form F(|D| , I) 2 for A at a bouncing ball orbit 7. The normal form lives on the model 
space fi := [0, L] s x K y , which carries a natural abelian algebra of pseudodifferential 
operators A = (\D S \,I), where: \D\ = \D S \ = y/—D*D s with Dirichlct boundary 
conditions on [0, L] and where / is a quantum action operator. Here, D x = \-§^- In 
the elliptic case, I = i(|Z?| _1 Z3^ + |£% 2 ) is the 'transverse' homogeneous harmonic 
oscillator, while in the hyperbolic case I = ^^D]' 1 D 2 — \D\y 2 ). By definition, the 
normal form of y/~K is a first order polyhomogeneous symbol 

(i) Fm t ) „w + °t + !^ + »M + ... 

in \D\, I with pj a polynomial of degree j + The coefficients B^ mentioned above 
are the coefficients of these polynomials (cf. Q | Z.1 | |Z.4|). 



Heuristically, v A should be microlocally conjugate to I) near correspond- 

ing bouncing ball orbits. This should imply that y/A and F(\D\,I) have the same 
wave invariants at iterates of their corresponding bouncing ball orbits, and conse- 
quently that the formal function F is a spectral invariant. As may be anticipated 
the boundary gives rise to many complications and we cannot quite implement this 



SPECTRAL DETERMINATION OF ANALYTIC BI-AXISYMMETRIC PLANE DOMAINS 3 



outline. In fact, we work as much as possible in the open space containing the 
bounded domain. The normal form is actually used to define a parametrix for the 
Dirichlet wave kernel near 7. 

1.1. Outline of the proof. The main contribution of this paper is the construc- 
tion of a normal form F(\D\,I) of A around a bouncing ball orbit 7, and the proof 
that the normal form is a spectral invariant. Here are the main steps. 

1.1.1. Straightening the domain. To conjugate to the normal form, we first use 
a special map <E> introduced by Lazutkin in JD to 'straighten the domain' near 
AB, i.e. to carry an open neighborhood of 7 in Q to an open neighborhood of 
[0, L] x {0} in f2 G . In fact, <& will be defined in a neighborhood in R 2 of AB. 
In the elliptic case, Lazutkin's map additionally puts the metric into the normal 
form ds 2 + b{s, y)[y 2 ds 2 + dy 2 }, or equivalently the Laplacian has the normal form 
D 2 + B(s, y)(y 2 D 2 + D 2 ) + LOT (lower order terms). We will modify Lazutkin's 
construction so that in the hyperbolic case the Laplacian has the normal form 
D 2 + B(s, y)(—y 2 D 2 + D 2 ) + LOT. We do not use the metric normal form in this 
paper beyond the quadratic term and we do not really need to define the straight- 
ening map by Lazutkin's method. We do so anyway because there is no advantage 
to constructing another map and because we believe the details of Lazutkin's con- 
struction could be useful in the general inverse problem. 

The straightening map carries A to a variable coefficient Laplacian A in a neigh- 
borhood of [0,L] x {0} in M 2 . A given wave invariant depends only on a certain 
germ of the Laplacian at the orbit, and consequently only on a certain germ of the 
straightening. So with no loss of generality we may assume that A is a polynomial 
differential operator in the y variable with analytic coefficients in s. We define the 
Dirichlet Laplacian Afi o to be this operator with Dirichlet boundary conditions on 
dn o . 

1.1.2. Conjugation to normal form. The next step is roughly to conjugate Aq o 
to a microlocal normal form near 7 D := [0, L] x {0} by an FIO (Fourier integral 
operator) on the model space. Intuitively, we would like to construct a microlocally 
invertible FIO W which 'preserves Dirichlet boundary conditions' and such that 
VF _1 AW / ~ F(|D|, I) 2 modulo an 'acceptable remainder' near 7 . 

Let us be more precise. First, by 'preserving Dirichlet boundary conditions' we 
mean that W should carry the domain of the Dirichlet Laplacian A as an unbounded 
operator on fl a to the domain of F(|£>|, I) 2 as an unbounded operator on f2 Q . Most 
significantly, Wu = on 9f2 G if u € i?g(Q ). 

Second, let us be more precise about the sense in which we are conjugating A to 
F{\D\, I). It is technically complicated to conjugate a boundary value problem by 
Fourier integral operator methods, so we do something simpler which is sufficient for 
the proof of our theorem. Namely we observe that both the Dirichlet wave operator 
cos^V An) and the normal form wave operator costF(\D\, I) are restrictions to their 
domains of well-defined Fourier integral operators in microlocal neighborhoods in 
the open space of the bouncing ball orbits . We can use the intertwining operator 
W on the open space to 'pull back' cos tF(\D\, I) (or more precisely its 'odd part') 
to a parametrix for cos ty/ An) in the interior of fl . It is only in this weak sense 
that we conjugate the wave group to normal form. Since the wave trace Tr cos t\fK 
at 7 involves the wave kernel only in the interior, we can compute it in terms of 
this parametrix and hence in terms of the normal form. 
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Our goal then is the construction of an intertwining operator to normal form 
in this weak sense and the characterization of the error term in the conjugation. 
There are two main issues we would like to emphasize in this introduction. The 
first has to do with solvability of the conjugation equations. Those familiar with 
the conjugation to normal form in the boundaryless case will recall that (just as in 
the classical conjugation to Birkhoff normal form) it is based on solving a sequence 
of homological equations [Q, \D\ + j-I] = KNOWN for the infinitesimal intertwining 
operator Q. This is a first order equation for the symbol of Q and it may seem 
mysterious that one can solve these equations with two boundary conditions (one 
at each boundary component). The second point to explain is the relevant notion 
of 'acceptable remainder'. 

1.1.3. Acceptable remainders. To clarify the second of these points, we recall that 
the link between the spectrum and normal form is through the coefficients in the 
singularity expansion 

(2) 

OO 

TrU(t) = Cy(t - L 1 + iOy 1 + a 7 o log(t - L 1 + iO) + a lk {t - L 1 + i0) k \og(t - L 7 

fe=i 

of the trace of the wave group U(t) — exp(itvA) at t = L 7 . When L 7 is the 
length of a bouncing ball orbit (or in general a periodic reflecting ray), the wave 
trace expansion is very similar to the boundaryless case in that the singularity is 
Lagrangean and the coefficients may be calculated by the stationary phase method 
(see Corollary (2.5)). In another language (cf. Corollary (2.6)), the wave invariants 
a 7 fc are non-commutative residues res{4^U (t)\ t= L^ of the wave group and its time 
derivatives. As already proved in [^) |Z.1|, it follows that only a certain amount 



of data from the Taylor expansion of the symbol of A along 7 goes into a given 
wave invariant 7. The precise statement is that a 7 fc depends only on the class of 



a A modulo the symbol class S 2 ' 2 ( k+ ^ (V, R7) of Boutet de Monvel pl\| . Here, V 
is a conic neighborhood of the symplectic cone R7. The bigrading of symbols is 
in terms of symbol order and order of vanishing along 7 (cf. §2). Terms of low 
symbolic order or of high vanishing order along 7 do not contribute to a 7 fe (cf. 



Proposition (2.7)). Thus we need to construct an FIO W which preserves Dirichlct 
boundary conditions and which conjugates Aq — > _F(|Z) S |, J) 2 modulo a remainder 

1.1.4. Reduction to a semiclassical problem. To do this, we convert the problem to a 



semiclassical conjugation problem as in [Z.l]. With the proper semiclassical scaling 
of symbols, elements of S 2,K are detected by their coefficient in the semiclassical 
parameter TV -1 . The goal then is to construct a semiclassical intertwining operator 
Wn preserving Dirichelt boundary conditions and conjugating the scaled version 
of A to normal form modulo a sufficiently high power of TV -1 . 

As in |ZT[ and elsewhere, we construct W N as a product W N = 11°^ e^' \P+iQ)i/* 



as a product of elliptic semiclassical pseudodifferential operators. The exponents 
will just be Weyl pseudodifferential operators (P + iQ)j/2(s,y,D y ) on the trans- 
verse space R; here, P, Q are assumed to have real-valued symbols. The real part 
P is of two lower orders in (y, D y ) than is Q. As mentioned above, the condition 
that Wn intertwines to normal form translates into homological equations for the 
exponents. The reason why we can solve these equations while preserving Dirichlet 
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boundary conditions is that the boundary condition only affects the 'odd terms' in 
Pj/2 and the 'even terms' in Qj/2 with respect to the involution (y,rf) — > (y,—rj). 
The equations for the even/odd parts of Pj/2, Qj/2 are coupled except for 'diagonal' 
terms which are powers of /. Hence one can eliminate the 'non-diagonal' even parts 
to get second order ordinary differential equations for the odd parts. The boundary 
problem for these can be solved in the hyperbolic case, and in the elliptic case they 
can be solved as long as a is independent of n over Q. The remaining powers of / 
constitute the terms in the normal form. 



1.1.5. Spectral invariance of the normal form. The intertwining operator W con- 
jugating An to F(\D\,I) 2 modulo g 2 > 2 (fc+ 2 ) can be used to construct a microfocal 
parametrix for the Dirichlet wave group modulo similar acceptable errors. The 
definition of the parametrix is £(t) — WF (t)ipvln W-i, where W-\ is a microfo- 
cal inverse to W, where F a (t) is the odd part of the normal form wave group 
costF(\D\, I) (the odd part satisfying Dirichlet boundary conditions at s = 0,L), 
where 1q o is the characteristic function of f2 and where ipv is a microfocal cutoff 
to a suitably small conic neighborhood of j . The wave invariants of the Dirich- 
let Laplacian Aq can then be calculated in terms of the normal form coefficients, 
and conversely ( as in Q (see also |p5.l[ ||Z.2f| ) the normal form coefficients can be 
determined from the wave invariants. Hence, the coefficients of the normal form 
F(\D\,I) are spectral invariants of An; a fortiori, the classical normal form of the 
Poincare map "P 7 is a spectral invariant. 



1.1.6. Conclusion of the proof. Theorem (1.1) then follows from a theorem of Colin 
de Verdiere |CV| that a bi-axisymmetric analytic plane domain is determined by 
the Birkhoff normal form of of a non-degenerate elliptic axial orbit 7. His proof 
works as well when 7 is hyperbolic and therefore we can conclude the proof in either 
case. 



1.2. Future problems. It should be remarked that F(\D\,I) is explicitly con- 
structed by the algorithm of this paper. In conjunction with Lazutkin's construction 
of a metric normal form, one can get explicit albeit complicated formulae for the 
normal form coefficients as polynomials in the Taylor coefficients of the boundary 
defining functions at the points A, B. In the future we plan to take up the obvi- 
ous question of whether one can determine the Taylor coefficients from the normal 
form coefficients. Colin de Verdiere's theorem shows that the principal symbol of 
the normal form alone is enough to determine these coefficients when the Taylor 
coefficients at A and B are the same and when the odd coefficients at A vanish. 
In less symmetric cases one will have to go into lower order terms in the normal 
form. In the simpler but somewhat analogous case of surfaces of revolution, it was 
necessary to go two steps below the principal symbol level to determine all of Taylor 
coefficients from the normal form. Rotational symmetry is analogous to left-right 
symmetry in a plane doman, so we suspect one can solve the inverse spectral prob- 
lem at least for left-right symmetric analytic domains using just the wave invariants 
at one orbit. 



1.3. Acknowledgements. The research in this paper was partially done dur- 
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in the hyperbolic case. 

2. Preliminaries 

In this section we establish some standard notation and terminology concerning 
the wave equation with mixed boundary conditions on a bounded smooth domain 



and its associated billiard flow. For further background we refer to |GM| [PS] 



2.1. Billiards on plane domains. First we recall the definition of a non-degenerate 
bouncing ball orbit for billiards on a smooth domain in W l 

Let C K n be a smooth domain, and let dfl denote its boundary. Let g be a 
smooth metric defined in a neighborhood of O. The billiard flow $* of (O, g) is the 
geodesic flow of g in the interior T*(intiY), with the usual law of reflection at the 
boundary (i.e. the tangential component of the velocity remains the same but the 
normal component changes its sign). Here, intCl denotes the interior of Q. 

From the symplectic point of view, there are two important hypersurfaces in 
T*(f2): the manifold Tg^fl of (co) vectors with footpoints on d£l and the unit cotan- 
gent bundle S*ft. We understand by T*VL the restriction T*M. 2 \q. The intersection 
Sgofl — TqqVI n S*fl is transversal. The characteristic foliation of S*tt is spanned 
by the generator of the geodesic flow. The projection ir : Tg n fl — > T*fl which 
projects a (co)vector at x £ <9£1 to the (co)tangent hyperplane defines a real line 
whose fiber at x is the normal bundle N x (dQ). We denote by v x the inward unit 
normal so that N x (dfl) = M.i/ X . This line bundle coincides with the characteristic 
foliation of Tg n f2. The image of S^fl under tt is the unit disc bundle D*(dfl). 
We will identify it with the space Sg^fl of inward unit tangent (co)vectors with 
footpoints on dfl. 

One then defines the billiard ball map (3 : D*(dfi) — ► D*(dfl) as follows: lift a 
tangent (co) vector v to d£l of length < 1 to SqqFI and move along the corresponding 
geodesic until the ball hits the boundary. Then project its tangent vector to D*dtt. 

By a reflecting ray one means a broken geodesic or billiard trajectory of the 
billiard flow whose intersections with the boundary are all transversal. We denote 

the successive points of contact with the boundary by qo, qi, q2, Of special 

importance here are the periodic reflecting rays where q n — qo for some n > 1. We 
will denote such a periodic trajectory by 7 : [O, T] — > S*Q. By a bouncing ball orbit 
7 one means a periodic reflecting ray where <?o — 92- The projection to Q consists 
of a segment <?o9i which is orthogonal to the boundary at both endpoints, i.e. go9i 
is an extremal diameter. The period is of course twice the length of the segment, 
which we denote by L. We write go = A, q\ = B. 

The Poincare map "P 7 of a periodic reflecting ray 7 is the first return map to a 
symplectic transversal. A natural symplectic transversal is given by a neighborhood 
of 7(0) in Sgftfl. Thus for v S Sg^fl, V-y(v) is obtained by following the broken 
geodesic thru v until it reflects from <9f2 the nth time in an inward vector. The 
linear Poincare map P 1 is defined to be tfP 7 (7(0)). 

Equivalently, one can identity Sgi-fil with D*{yi) as above and define V 7 as a 
map on a neighborhood S 7 of the zero vector at an endpoint of 7 in D*(il). Then 
for a periodic reflecting ray 7 with n reflections, P 7 may be identified with /3™|s . 

Definition 2.1. A bouncing ball orbit 7 is said to be non- degenerate if the eigen- 
values of P 1 are not roots of unity. There are two cases: 7 is 
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(i) non-degenerate elliptic if the eigenvalues of P 7 are of the form {e ±ta } with 
a/ it i Q.; 

(ii) non-degenerate hyperbolic if the eigenvalues of P 7 are of the form {e ±A } for 
some A G K+. 



For background on bouncing-ball orbits we refer to [B.B] (see §5.2). If we denote 
by Ra, resp. Rb the radii of curvature of the boundary at the endpoints A, resp. B 
of an extremal diameter, then one finds that ellipticity is equivalent to the condition 
that L < Ra + Rb,L > Ra, L > Rb or else that L < Ra + Rb,L < Ra, L < Rb- 
If AB is a local minimum diameter than L < Ra + Rb while if it is a local 
maximum diameter then L < Ra + Rb- So a (non-degenerate) local maximum 
diameter must be hyperbolic; a local minimum diameter is elliptic if it satisfies 
the additional inequalities above. Under our symmetry assumption Ra — Rb, so 
ellipticity is equivalent to the statement that the center of curvature at A lies below 
the horizontal axis. 

2.2. Wave trace on a manifold with boundary. Let A denote the Laplacian of 
a metric g defined in an open neighborhood of a bounded smooth domain Oct 2 . 
The Dirichlet Laplacian Aq of fl is then defined to be the self-adjoint operator on 
L 2 (fl) with domain 

(3) Vom(An) = {ue H^(n) : Au £ L 2 (n)} 
where i?Q (SI) denotes the closure of Co°(fi) under the norm 

h«ii;=£ii£i£" 

j<2 3 

and where where Au is taken in the sense of distributions. We recall that Hq(CI) = 
{u G ff^R 2 ) : supp u C ft}. 

Now let E(t, x, y) be the fundamental solution of the mixed wave equation with 
Dirichlet boundary conditions: 

^ = AE onMxOxO 

(4) E(0,x,y) = S(x-y) |f(0,z,y) = 
E(t,x,y)=0 (t, x, y) G M x dfl x Q. 



As discussed in ( [|GM| , §5]|PS|), the fundamental solution is not globally a La- 
grangean distribution. However, for any T > there exists a conic neighborhood 
Fx of the bouncing ball orbit 7 so that the microlocalization of E(t, x, y) to Tt is 
Lagrangean for \t\ <T. More precisely, 



Theorem 2.2. (]GM| ; Theorem 4-1 and Proposition; or [PS], $6) Let (x,£) € 
T*(intfl) or let x G dfl and suppose that £ is a non-glancing (co-) direction at x. 
Then there exists a conic neighbhorhood O of (x, £) in T*R™ — and a Fourier 
integral distribution V(T, x, y) essentially supported in O such that for any pseu- 
dodifferential operator \o{x,D) essentially supported in O, xo{x, D)V (T,x,y) — 
Xo(x,D)E(T,x,y) <E C°°(R x ft x O). 

By a partition of unity, this leads to a Fourier integral formula for the wave trace 
near a periodic reflecting ray: 
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Theorem 2.3. [ A.M] Suppose that T is an isolated point of the length spectrum of 
with the following properties: 

(i) T^\dQ\; _ 

(ii) All of the closed billiard trajectories of length T are non- degenerate reflecting 
rays. 

Then modulo smooth functions in t near t — T we have: 
Ea^gSpCVa) cosA ^ = In ^C*' x > x ) dx 

= E J I a V^ ) (t,x,x)dx 

where V± (t, x, y) are Fourier integral distributions associated to the broken geodesic 
flow near the closed trajectories of length T . 

The wave invariants can therefore be obtained by applying the method of sta- 
tionary phase to an oscillatory integral. An apparent obstruction is that the domain 
of integration is a manifold with boundary. However, Guillemin- Melrose prove: 

Lemma 2.4. Let il C W l be a smooth domain with boundary and let <j) E C°°(R") 
have clean critical point sets. For each critical value X let C\ denote the critical 
set in K™. Suppose that for each X, C\ intersects dfl transver sally. Let U\ be a 
neighborhood of C\ with the properties: 

(i) U\ contains no critical points of (f> except C'\; 

(ii) U\ n dQ, contains no critical points of 4>\dn except C\ n dfl. Then the usual 
stationary phase expansion is valid: 

oo 

a(x)e lT ^ x) dx ~ e lTX T~ k/2 V a,r'\ (r -» oo) 

where k is the codimension of C\ in R™ and where on are integrals of a and its 
derivatives over C\ D Q. 

2.3. Wave trace invariants as non-commutative residues. We may summa- 
rize the relevant result on the Poisson formula as follows: 

Corollary 2.5. Let 7 be a non- degenerate billiard trajectory whose length L 7 is 
isolated and of multiplicity one in Lsp(Q). Let Tl be a sufficiently small conic 
neigbhorhood o/]R + 7 and let %p be a microlocal cutoff to Tl- Then for t near L 7 , 
the trace of the wave group has the singularity expansion 

00 

Tr^pE{t) = c^(t-L 1 +iO)~ 1 +a 7 o log(t-£ 7 +iO)+V] a jk (t-L 1 +iO) k log(£-L 7 -H0) 

fc=i 

where the coefficients a 7 fe are calculated by the stationary phase method from a 
Lagrangean parametrix. 

This corollary allows us to identify the wave invariants as non-commutative 
residues as in |G||Z.l]. Recall that if A is a Fourier integral operator, and if P 



is any positive elliptic first order pseuodiffcrcntial operator, then the zeta function 
C(z, A, P) := TrAP~ z has meromorphic extension to C with at most simple poles. 
The residue at z = is referred to as the non-commutative residue res(A) of A. 
It is independent of A and is a tracial invariant, i.e. res(WAW~ 1 ) — res A. In 
the boundaryless case, one has a 7 fe = res{j^) k E(t)\ t= L^ ■ The only ingredients in 
the proof are the Lagrangean property of E and a canonical transform between 
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TrAe itp and TrAP~ s . Hence the same result remains valid in the case of periodic 
reflecting rays of the boundary case: 

Corollary 2.6. If '7 is a periodic reflecting ray, then a lk = res( j^) k E(t)\ t= L i . 



Now we recall some elementary results from (g| and ZJ] on the data of the 



domain or metric which go into a given wave invariant a 7 fe in the boundaryless 
case. Analogous results hold in the boundary case, but we postpone stating them. 

In the following, let P be any first order pseudodifferential operator of real 
principal type on a boundaryless manifold U and assume for simplicity that all 
closed orbits of its bicharacteristic flow are non-degenerate. Let a 7 fc(P) denote the 
kth wave invariant at a closed orbit 7 for e ltp . Since it is calculated by a stationary 
phase expansion at 7, it is obvious that a 7 fc(P) depends only on certain part of the 
jet of the complete symbol of P around 7. 

To state the precise result, let us Taylor expand each term in the complete (Weyl) 
symbol p(s, er, y, ry) ~ p± + p a + . . . of P at M. + j: 

Pj{s,v,y,v) = ct j Pj(s, -) = o- J (p [ ° ] +pf ] + ...) 

with p^ m ' (s, 1, y, the part which is homogeneous of degree m in (y, Set Pj :— 
Opipj), Pj m] := Op{pf ] ) and Pf N = ^ m ^ N P\ m \ Then we have: 



Proposition 2.7. (Jzj|, Proposition 4.2) a lk {P) = a 7fc (Pp (fc+2) + Pr 2(fc+1) + 
■■■P°-k-i)- 

Thus, it is sufficient to define the normal form and the intertwining operator in 
a microlocal (conic) neighborhood V of R + 7. In (s,a,y,ri) coordinates, we may 
define the cone by: 

(5) V = {(s, a, y, n) G T*U e : | V | < e, \ V \ < e\a\}. 

Let ipv be a microlocal cutoff to V , with symbol identically equal to one in a slightly 
smaller open cone around M7 and put A e :— ip e AQip e . Ay and A of course have 
the same microlocal normal form around 7 so for notational simplicity we often 
drop the subscript and leave it the reader to recall that the operator is cutoff. We 
will re-instate the cutoff at the crucial point of calculating the residue. 

Since terms which vanish to too high order at 7 or which have too low a pseu- 
dodifferential order do not contribute to a 7 fc we introduce a bi-grading on symbols 
in terms of order as a symbol and order of vanishing along 7. Let R + 7 C T*U 
be the cone thru an embedded curve 7 € T*U and OjS m (U) denote th e cla ss of 
symbols of order m over U which vanish to order j along 7. Following (|BM[|, see 



also [BMGH]) we denote by S m ' k (V, 1R+7) the the class of symbols microsupported 



in V which admit asymptotic expansions 

00 

(6) a ~)^a m -j, a m -j{x,r£) = r m_:, a(a;, <^), a m -j G O k ~2 3 S m ~ J . 
j=o 

Here, k G N and is no condition if 2j > k. We denote by Op w S m ' k the Weyl 
pseudodifferential operators with complete symbols in S m,k . 
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In local coordinates (s, y) with 7 = {y — 0} and with dual symplectic coordinates 
(cr, 77), a e S m ' k (V,M+j) if 

a{s,a,y,rj) ~ J2%o a m-j,fe(s, cr, y, 77) - ^°^ o m -'-> a m ^ oM (s, 1, y, 77/0-) 

(7) 

a m -j >k (s,l,y,r)/a) S O k - 23 S°. 

The coefhcient a m -j(s,l,y,r)/cr) is homogeneous of degree and the assumption 
that a is microsupported in becomes that a m „j is supported in the transverse 
baU£ e = {|(i/W)|<e.} 

Any symbol a £ S m (V) may be expanded as a sum of symbols in S m,k . Indeed, 
let a m -j-k be the term of degree k — 2j in its Taylor expansion for k > 2j 

(8) a m _-, ; fe = ^ am-j;k-2j,a,p{s, )y° {v / <J )' 3 ■ 

\a\ + \p\=k-2j 

Then the asymptotic sum a^ rnk ^ ~ Y/j—^ a m—y,k belongs to S m ' k (sharp) and 
a = Sfc°=o a (m,k) ■ I n particular, we may expand A in this form, and a lk depends 
only on the class of A modulo S 2,k . 

We summarize the discussion by restating Proposition ( |2.7| ) in terms of these 
symbol classes: 

Corollary 2.8. Suppose that A, B E Op w S 1 {V) and that A = B mod OpS x ' k {V). 
Then a 7 fc(A) = a^fc(S). 

3. Straightening the domain 

As discussed in the introduction, the normal form lives on the model domain 
fl = [0, L] x R, or more precisely in a microlocal neighborhood of T* [0, L] in its 
cotangent bundle. In this section we introduce the analytic objects on the model 
domain and explain how to transfer A, in a neighborhood in R 2 of a non-degenerate 
(elliptic or hyperbolic) bouncing ball orbit 7, to a variable coefficient Laplacian A 
in a neighborhood of [0, L] x {0} in the model. We emphasize that all maps and 
operators that we discuss in this section extend to open domains containing the 
various manifolds with boundary. 

3.1. The model domain fl Q . The configuration space of the model is the infinite 
strip Q a . We denote the coordinate on [0, L] by s and that on K by y , with dual 
cotangent coordinates cr, 17 on T*[0,L] x T*R. We also denote by I e — \{t} 2 + y 2 ) 
the elliptic action variable and by Ih = ^(r] 2 ~U 2 ) the hyperbolic action variable on 
T*K. To simplify notation we often just write I for the relevant action variable. In 
the Poisson algebra of T*(K 2 ) we consider the maximal abelian subalgebra A c i — 
(\a\,I). The model (classical) Hamiltonians are those of the form H a = \a\ + j^I 
which generate linear Hamiltonian flows. As in the case of straightline motion, H a 
generates a broken Hamiltonian flow on T*[0,L] x T*M. when equipped with the 
boundary condition that the trajectory is reflected by r when it hits the boundary. 

We will view tt as a submanifold with boundary of the 'open space' S\ L x R 
with S\ L = R/2LZ = [—L,L]. Many of our operators will be supported in a 
neighborhood U e — (— e, L + e) x M.) of f2 D in S\ L x R. We denote by r e (s) a smooth 
cutoff to U e with t c = 1 in a smaller neighborhood of f2 . 
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Since o , U e , and S\ L x R are products, their algebras of pseudodifferentials 
operators are easily described in terms of pseudodifferential operators along S\ L 
and transverse pseudodifferential operators on R. 

In the direction of S\ L , we introduce the algebra ^*{S\ L ) — (s, D s ) of standard 
pseudodifferential operators on S\ L . We distinguish the element \D a \ = y/—D* s D s 

■ Tvks 

with eigenfunctions e % ~ associated eigenvalues \k\. We also define the subspaces 
and projections: 

H 2 ±(Sh) = ©r = oCe ±l ^ n± : L 2 — > H% 

n +fc : L 2 {S\ L ) -> Ce^ II_ fc : L 2 (S\ L ) -> Ce"^ 

(9) 

io^2L) = ©r=oCsin(^M) n o :L 2 ^L^ n o /( S ) = - r/( S )) 

Here, r/(s) = /(— s) where — s is taken modulo 2L, i.e. r is reflection through the 
boundary. We may tensor the subspaces with L 2 (R) to get corresponding subspaces 
of and operators on L 2 (S\ L x R) and we use the same notation for these. 

We will sometimes identify functions on Q with odd functions on S\ L x R. More 
precisely, let us put: 



(10) 

n fe :L 2 ([0,L] xR)^Csin(42)(g)L 2 (R) 

l [oM (s):L 2 (Sl L )^L 2 ([0,L}) 

A:L 2 {[0,L])^L 2 odd {Sl L ) 



l[o,L]( s ) — the characteristic function of [0, L] 
Af = the odd extension of / 



A* ■■ L 2 odd {Sl L ) -+ L 2 ([0, L]) A*g = l [0 . L] n o g 

\D s \:H^[0,L])^L 2 ([0,L}) \D.\ sin(^M) = f sin(2M) 

We use the notation \D S \ (or simply \D\) simeltaneously for the 'Laplacian' on 
S\ L and the Dirichlet Laplacian on [0, L}. No confusion should result since they are 
defined on different domains, and moreover the definitions are compatible under 
the above identification, as the following proposition shows. 

Proposition 3.1. We have: 

(i) A : L 2 ([0, L]) - L 2 odd {Sl L ) and AA* = U , A* A = Id. 

(ii) A*\D S \A = \D S \ 

The proof is obvious so we omit it. 

In the transverse direction, we first introduce the isotropic Weyl algebra W* . 
This is a completion of the algebra £ :—< y,D y > of polynomial differential op- 
erators on R. We denote by £ n denote the subspace of polynomial differential 
operators of degree n in the variables y, D y . We also denote by £™ the polynomials 
all of whose terms have the same parity as n. In the isotropic Weyl algebra W*, 
the operators y, D y are given the order i, so that 

(11) £ n cW n/2 , [£ m 1 £ n ]d£ m+n - 2 . 
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We note that such operators are not standard homogeneous pseudodifferential op- 
erators, but can be rescaled to this form. The rescaled algebra is generated by 
the first order homogeneous pseudodifferential operators 2/|D s |2 and \D S \~^ D y . In 
the open space the relevant algebra is the algebra of homogeneous pseudodifferen- 
tial operators over U e generated by D s ,y\D s \2 , \D S \~? D y . To be more precise, this 
construction is only well-defined in a microlocal neighborhood of 7 where |D S | is 
elliptic. 

It should also be recalled that a w (y,D y ) g £ n with real- valued symbols are es- 
sentially self-adjoint operators on S(J3L) (the Schwartz space), i.e. have a unique 
self-adjoint extension to L 2 (M.). Hence their exponentials e la ™^ v ' Dv ^ are unambigu- 
ously defined. 

3.2. Half-density Laplacian. In order to deal with self-adjoint operators with 
respect to the Lebesgue density dsdy, we pass from the scalar Laplacian to the 
(unitarily equivalent) 1/2-dcnsity Laplacian 

Ai := JV2AJ-1/2 

2 

(12) =El =1 J- 1/2 D x ^JD X] J-^ 

= g^Dl + g 22 D 2 y + 2g^D s D y + l^D s + \T 2 D y + a . 

Here we write (s, y) — {x\, X2), D Xj — j§^-- The functions Tj are real valued. Since 
it is self-adjoint relative to the Lebesgue density, its complete Weyl symbol is real 
valued. Henceforth we denote the I/2-density Laplacian simply by A. 

3.3. Straightening the domain. Let us now explain how to transfer the Lapla- 
cian to the model domain. As in the boundaryless case, the model space is in some 
sense the normal bundle of the orbit. This is literally correct in the boundaryless 
case and the exponential map along the normal bundle N 7 of 7 can be used to trans- 
fer A to the normal bundle. In the case of a bouncing ball orbit in the boundary 
case, the normal bundle and exponential map are ill-defined at the reflection points 
but Lazutkin has constructed a nice replacement for them. Namely, he constructs 
a map $ which straightens the domain to a strip near a stable elliptic bouncing ball 
orbit and which simeltaneously puts the Laplacian into a preliminary normal form. 
We will modify his method to encompass hyperbolic bouncing ball orbits as well. 
As mentioned above, we do not need the full details of the map or metric normal 
form here. Hence we only sketch the construction of $, referring the reader to (l| 
for the details. In the following, f2 e denotes an e-neighborhood of AB in il. 

Definition 3.2. By a transversal power series map from £l e to tt a we mean a 
formal power series 

$:fi e ^[/ £! $(s,y) = (s,y) 

of the form 

S = S + J2m=2 K mm(s)y m 

y = T,%Lii> P p(s)y p 

with real valued analytic coefficients in s extending analytically to a neighborhood 
of [0, L] in C and satisfying the boundary conditions: 

{s = 0} U {s = L} = $(dfi e ). 
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We use the language of formal power series since we only need to use a polynomial 
part of the map to construct the normal form up to a desired accuracy. Indeed, a 
given wave invariant a 1 k only involves the 2k + 4 -jet of <&. The convergence of the 
series is irrelevant to our purposes and we will not discuss it. The following Lemma 
was in effect proved by Lazutkin in jE) in the elliptic case. 

Lemma 3.3. Suppose that AB is a bouncing ball orbit. Then there exists a transver- 
sal power series map $ : (r2 e ,<90 £ ) — ► (Q ,dQ ) which straightens the domain and 
puts A in the form: 

A ~ D| + B(s, y)(y 2 D 2 + D?) + T' S D S + T' y D y elliptic case 

A ~ D\ + Bis, y){D 2 y - y 2 D 2 s ) + r' s L> s + T' y D y hyperbolic case 

in the sense that the left and right sides agree to infinite order at y = 0. Here, 
B(s,y) is a transversal power series. 

Of course, we only need this form of A to order K construct the normal form 
modulo OpS 2,K , and only in the principal terms do we need to know the exact 
form. Therefore we only briefly recall the proof of the lemma and only discuss the 
principal terms in detail. 

Under any map the usual (scalar) A conjugates to 

(13) 

A := $*-!A$* = \Vs\ 2 D 2 + \Vy\ 2 D 2 + 2(Vs,Vy)D s D y + -AsD s + -AyD y . 

" ii 

In the case of a transversal power series map, the coefficients are also transversal 
power series of the form: 

\Vs\ 2 = l + E™=2am m (s)y m 

|Vy| 2 = E™=oWs)y m 
(14) As = E" =0 Cm m (s)y m 

^y = HZ= d m m{s)y m 

(Vs,Vy) = Em=l e mmO)y m 

The volume density in the new coordinates, J(s,y) has a similar form. The 1/2- 
density Laplacian ($* _1 A$*)j. in the transformed coordinates is given by 

(15) 

A := J-^($*- x A$*) J? = \Vs\ 2 D 2 + \Vy\ 2 D 2 + 2(Vs, Vy)D s D y + -T S D S + -T y D y + K 

a ii 

where 

T g = -| Vs| 2 J. log J - (Vs, Vy)^ log J + As 
(16) Ty = -|Vy| 2 ^logJ~ (Vs,Vy}JUog J + Ay 

K = J-2($*- x A$*)(J5). 

Combining with the power series expressions in ( p^ ) one has an expression for A 
modulo terms vanishing to order K at y = 0. 
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The principal term (modulo S' 2 ' 3 ) is given by D 2 + (a22(s)y 2 D 2 + en(s)yDy + 
boo(s)D^. We now choose the coefficients in $ to put it in the preliminary normal 
form of 

D 2 + \b 00 {s)[y 2 Dl + £>?] elliptic case 

(17) 

D 2 + |6oo(s)[-D| - y 2 D 2 } hyperbolic case 
This requires the coefficients to solve the equations: 

• en = 

(18) • a,22 = boo elliptic case 

• 0,22 — —boo hyperbolic case. 

To solve these equations we use the following expressions from (jy, (4.5) - (4.9)): 

022 (s) = (2k' 22 + Ak1 2 )^ 2 

(19) M*)=i&(«) 

coo(s) = 2k 2 2(s) 

Since they are obtained purely algebraically, the same equations hold in both the 
elliptic and hyperbolic cases. They are easiest to solve if we make the substitution 
0n — ipii ■ One then has: 

eu =0 = 

(20) a.22 = b o => O'd = p- elliptic case 

022 = —boo Q'u = — W hyperbolic case 
The solutions have the form: 



^ f ^ 2 elliptic case 

(21) ' 

On = \J-l + hyperbolic case 

for some constants I, sq. To verify this it is easiest to substitute £ = 6* n . The 



equation for £ is then h^" = t(£') 2 ± 1 (+ = elliptic, - = hyperbolic). Solving for 



&oo gives 



b Q0 (s) = ^ { l- S0 )2 (elliptic), b 00 (s) = - - * g _ So)2 (hyperbolic) 

hence 
(22) 

j g ^ P ^ g . g 

boo(s)ds = - tan _1 ( — ^) (elliptic), y b o(s)ds = — tanh~ ( — ^) (hyperbolic). 
In the elliptic case one finds that 

(23) t = ^x {R A -s ) = y/(L-s )(R B -L + s ). 
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These conditions uniquely determine x a ,£ (see the pictures on p. 135 of [L.l]). It 

S-SQ 

also follows that coo(s) = — t 11 = r — zr. 

In §4.2 we will put the leading term into a canonical normal form with coefficients 
independent of s. 

3.4. Wave invariants revisited. Having straightened the Laplacian, and hence 
having transferred the information about the boundary into the metric, we can 
now state precisely just how much data of the boundary goes into a given wave 
invariant a 7 fc. Since the wave invariants at {y — 0} of A are calculated by the 
stationary phase method, we have, as in the boundaryless case: 

Proposition 3.4. a 7fc (A) = a 7fc (A^ 2(fc+2) + Af 2(fe+1) +■■■ + A° k _ 1 ). 

As A is a well-defined partial differential operator in a neighborhood of Q we can 
again reformulate the conclusion in terms of the symbol classes S 2 ' k (V, R + 7o) where 
7o denotes the bouncing ball orbit [0, L\ x {0}. The discussion in the boundaryless 
case remains valid, althought the microfocal neighborhood V now acquires two 
components: Since |oj > \ > in V, V = V+ U V- where V + = V n {a > 0}, V_ = 
V H {a < 0}. The two components are obviously interchanged by the canonical 
involution r(s, a, y, rj) = (s, —a, y, —r)) of T*([0, L] x R). 

Corollary 3.5. a 7 fc(A) depends only on the class of A modulo g 2 > 2 ( fc + 2 ). 

In view of the metric normal form, the data which goes in to the k-jet of A along 
{y = 0} is precisely the k + 2-jet of the function B(s, y) and hence the k + 2-jet of 
the boundary defining functions. Hence a 7 fc depends only on the 2k + 4-jet of the 
boundary at {y = 0}. 

We will actually need a slight generalization of ([3.5|) which is prove in precisely 
the same way. 

Proposition 3.6. Suppose that A G Op(S 2 ' 2 ( k+2) ) . Then resFAG = for any 
bounded Fourier integral operators F, G. 

4. SEMICLASSICAL NORMAL FORM 

As mentioned in the introduction, our approach is to convert the conjugation to 
normal form to a semiclassical problem. 

4.1. Semiclassical scaling. To introduce the semiclassical parameter, we make a 
semiclassical scaling of operators on the model space. Roughly, the scaling weights 
the tangential derivative D s by TV 2 , the normal derivative Dy by N and y by iV _1 . 
We define operators Tat, Mjv on the model space L 2 (il ) by 

. T N f(s,y):=Nf(s,Ny) 

(24) 

. M N f(s,y) :=e m s f(s,y) 

We then have: 

T* N DyT N = NDy 
(25) T* N yT N = N^y 



M* N D S M N = (N 2 + D s ) 
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Definition 4.1. The reseating of an operator a w (s,D s ,y,Dy) is given by 
a w N {s,D- s ,y,Dy) := M* N T* N a w (s, D s ,y, D V )T N M N . 

We have: 

Proposition 4.2. The complete symbol a% (s ', a , y , fj) of M^T^a w (s, D s ,y, D^TnMn 
is given by 

a%(s,a 7 y 7 fj)=a w (s,a + N 2 ,^y,Nfj). 
Proof The operator kernel of a w (s, D s , y, Dy) is equal to 

f f »(J^sl^(9-»')>* i((i " ,J '' )+(Fi '' >Wf '' 

Conjugating with Mjv amounts to adding N 2 (s — s') to the phase and hence to 
a translation a — > a + N 2 . Conjugating with T/v amounts to changing the phase 
(y — y 1 , fj) to jf{y — y', f]) and the amplitude to a(|(s — s'), a, — fj). Change 
variables fj — > iVry to get a Weyl pseudodifferential operator with amplitude a(^(s — 
s -'),a,i(f -|),^). " □ 

Definition 4.3. The semiclassically scaled (1/2- density) Laplacian is the Weyl 
pseudodifferential operator on Q defined by 

A N = M* N T* N AT N M N . 

In the straightened form, we have 

A N ~ (Ds + N 2 ) 2 + B(s, N-'y^fN^iD, + N 2 ) 2 + N 2 D 2 } 

(26) 

+{r g } JV [U 5 + N}+ N{T y } N Dy + {K} N 

where {f(s, y)}jv = T^f(s,y) = f(s,N~ 1 y) and where (as usual) the alternative 
± refers to the elliptic/hyperbolic dichotomy. 
The Weyl symbol of A has the simple form: 

(27) al(s,a,y,fj) := a 2 + B(s,y)I + K, K = A-1 

The linear terms vanish because they give the subprincipal symbol in the Weyl 
calculus and that of A equals zero. By the above proposition, we get upon rescaling 

(28) af N ~(a + N) 2 + B(s, ^y){±y 2 N~ 2 [a + N 2 ) 2 + N 2 V 2 ) + K(s, ly). 

4.2. Linearized problem. To get a sense of what is involved in putting A into 
normal form, let us first consider the 'linearized' problem which involves only the 
highest powers of TV: 

(29) al N =N A + 2N 2 [D S + b 00 (s)I] mod N 

where I denotes the quantum action operator: I e — \{D^ +y 2 ) in the elliptic case 
and I h = \{Dy^ — y 2 ) in the hyperbolic case. 

As in the boundaryless case [Z.l, §1] we can complete the conjugation of the 
linear/quadratic term to normal form by a moving metaplectic conjugation. When 
results apply mutatis-mutandi to both elliptic and hyperbolic cases, we denote the 
action simply by I. In the following proposition, ji denotes the metaplectic repre- 
sentation of SL(2, K) (strictly speaking, of its double cover but signs are irrelevant 
here). 
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Proposition 4.4. There exists an SL (2,R) -valued function a a (s) so that 

n(a a )*[D s + b 00 (s)I]iJ,{a a ) = TZ 

fj,(a a )(0) = n(a a )(L) = Id 
where TZ = D s + —I and where a = boQ(s)ds. 

Proof We first construct a function a(s) so that n(a)*[Dg + boo(s)I]iJ,(a) = Dg. The 
desired metaplectic operator is obviously given by 



(30) fi(a{s)) = exp(-i[ / b 00 (s)ds]I). 

Jo 

Hence a(s) = exp(— J/2 f£ boo(s)ds) where 

J = ( ® ^ q ^ elliptic case 

J=(J ° 1 ) hyperbolic case. 

The boundary values of a are given by a(0) = /, a(L) = cxp(— J J Q boo(s)ds) = 
exp(— al). Now let r a (s) = exp(^ Js) and put a a (s) = a(s)r a (s)^ 1 . Then /j,(a a (0)) = 
Id and fi(a)a(L) = exp(ai) o exp(— al) = Id. Further, /i(r Q (s)) conjugates D 8 to 
D s + f I. □ 

4.3. Semiclassical pseudodifferential operators. Semiclassical scaling produces 
a partial differential operator A n depending on a small parameter i /N and in the 
conjugation to normal form we will introduce other such operators. Let us pause 
to clarify the kinds of semiclassical pseudodifferential operators which will be of 



concern to us. Our discussion is based on ideas and notation from BMGH |. 

First, let us recall that the usual (admissible) semiclassical h — l/7V-pscudodiffcrcntial 
operators of order m on R™ are the Weyl quantizations 

a w (x, ±D; N)u(x) = (2irh)- n J J e^^a^x + y), N)u(y)dydt 

of amplitudes a belonging to the space 5™(T*M n ) of asymptotic sums 

oo 

a(x,t,N)~N m J2 a j( x >0 N ~ j 

3=0 

?2n\ 



with a,j € C°°(i? : 

Semiclassical scaling gives rise to symbols of the form a(s, a + N 2 , -^y, Nfj) where 
a is a polyhomogeneous symbol a ~ S^=o a m-j °f order m. As in (0) we may write: 



(31) a(s, a + N\ ifr Nfj) ~ £ N 2m ~ 2 i a m ^{s, I + jL, ±y, ±fj). 

j=o 

It is evident that semiclassical scaling produces symbols which behave in both 
transverse variables (y, fj) like standard semiclassical symbols in the fiber variable 
r\. Thus, scaled symbols are isotropic analogues of semiclassical admissible symbols 
in the transverse M. n . They also have a tangential dependence in 1 + ■£> which seems 
to have no precise analogue for standard semiclassical symbols. 
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Let i/j(s,a,y,r]) be the homogeneous cutoff to the cone V. Under rescaling it 

JV 7 ' N ' N ' 



goes over to the symbol N 2 ip(s, 1 + j?). For symbols independent of a or 



for a < eN the scaled symbol is supported in the transverse ball 

(32) V N :={(s,a,y,rj):\(±y,±fj)\<e}. 

It is useful to reformulate the condition that A <E OpS m,k in terms of the s caled 
symbol. Given A £ OpS m we define (with some modifications to [ BMGH ]) the 
formal differential operator 

(33) oft^A) = E ^(^) a (|)%(S.^».OS a ^ (|a|+l '' l+SW 
and put 

(34) °\ M M)= £ ^(^H^C 8 ^.*.^ 13 ?- 

\ a \ + \p\+2 3 =k 

Then a G 5 m ' fe (V r , R + j) if and only if is divisible by N~ k for (a;,f) = (0,0); 

hence (A) = a k - M) {A)N~ k mod A^^ 1 ) for (s,a,x,Z) G M+ 7o . We thus 

have 

Proposition 4.5. If a £ S m (V) then a £ S m ' k (V,R + j ) if and only if the formal 
Taylor expansion of a(s,N 2 , jfy,Nff) along K + 7 is divisible by N~( k+1 \ 

4.4. Conjugation to a semiclassical normal form. We now come to the princi- 
pal step in the conjugation to normal form: the conjugation of Am to a semiclassical 
normal form. 

We first conjugate by n(a a ) as in the linearization step to get the somewhat 
simpler form: 

(35) TZ N := n(a a )*A N n(a a ). 

We then wish to conjugate TZn to the semiclassical normal form 

(36) F N {lf ~N 4 + N 2 ^ + Pl (I) + N- 2 p 2 (I) + ■■■ 

by means of a semiclassical pseudodiffcrential intertwining operator Wn(s, y, Dy) 
which preserves Dirichlet boundary conditions. Thus the full intertwining opera- 
tor is /i(a a )WN- To avoid encumbring the notation we will also denote this full 
intertwining operator in (§5)by Wn- 

Let us now explain what we mean by semiclassical conjugation to normal form. 
We assume for simplicity that the bouncing ball orbit is elliptic, but the same argu- 
ment and result hold in the hyperbolic case. Roughly speaking, our object is to pro- 
duce bounded semiclassical pseudodifferential operators W k (s, y, Dy), W k (s, y, Dy) 
defined in the open neighborhood (even U e of Q and satisfying the following as- 
ymptotic relations on this domain: 

(i) AW+e™ k iD q (N k y) ~ F(k, q + \) 2 e™ k iW+ D q {N k y) 

(37) («) AW k -e-™ k iD q (N k y)~F(k,q+±) 2 er™ k iW k -D q (N k y) 
(Hi) W+e™ k iD q (N k y) - W~ e~™ k i D q {N k y) = G at s = 0, s = L. 
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Here, F(k,q + 5) = FN k {q + 5) with N k = \fk and D q is the qth normalized 
Hermite function. The precise meaning of ~ will be clarified below. The Laplacian 
A is the Laplacian acting on the open space U e . The condition (hi) implies that 
W+e™ k iD q (N k y) - e~™ k i D q (N k y) E H^(Q) so that it lies in the domain of 

Since CAC = A (with C the operator of complex conjugation), it suffices to 
construct W£ and to put = CW^C. Then, 

(38) W+e™ k iD q (N k y) - W~ e~™ k i D q (N k y) = 2&W+e™ k i D q (N k y) 
In order that 3H>+ e i7Tk i D q (N k y) lie Vom(A n ) it is thus sufficient that 

(39) CW£(0)C = W+(0), CW+(L)C = W+(L). 

This boundary condition (^) is correct in both the elliptic and hyperbolic cases, 
although the quasimodc construction which motivates it only works in the elliptic 
case. 

As a further preliminary, let us rewrite the equations ( |37| ) and the boundary 
conditions ([39]) in terms of the semiclassically scaled Laplacian. We observe that 
D q (N k y) = T Nk D q (y). Let us define: 

(40) W+=T£W+T Nk . 



Then (37)-(B9J) is equivalent to 



(i) A Nk W+D q (y) ~ F(k, q + \fW+D q (y) 
(41) (0) A Nk W^D q (y) ~ F(k, q + ±) 2 W^D q (y) 

(in) CW+(p)C = W+(0), CW+(L)C = W+(L). 

The operators = (s,y, Dy) will be essentially a family of pseudodiffer- 
ential operators on the transverse space, parametrized by s. Hence there are no 
subtleties involving the definition of pseudodifferential operators on manifolds with 
boundary. Moreover, is essentially applied only to a function of y. To be more 
precise, we first conjugate by /i(a a ) to put the linear term in normal form, and 
then is applied to a function in the kernel of 7Z. Hence we only require that 
the conjugation identity hold as operators applied to functions in the kernel of 7Z. 
Therefore we introduce the following notation: Given an operator A(s, D s ,y, D y ), 
we denote by A\ a the restriction of A to functions of y only, i.e. | G denotes the 
restriction to functions in the kernel of 1Z. 

The following lemma proves the existence of such a conjugating operator. We 
emphasize that the conjugation of A takes place over U e . 

Lemma 4.6. Let 1Z~^ = n(a a )* A.Nn(a a ) with — <^ Q in the elliptic case. Then 
there exist polynomial differential operators Pjj 2 (s, y, Dy) and Q™, 2 (s, y, Dy) of de- 
gree 2j + 2 on L 2 (R y ) with smooth coefficients in s £ (— e, L + e) and polynomials 
fj (I) of degree j + 2 such that: 

(a) Pj/2(s,y,r]) and Qj/2 (s, y, 77) are real-valued; 

(b) For each s, there is a formal N -expansion: 



DC 



(W+§ y ln N W Nf ~ ~ N4 + 2N2U + E N-in*_ l (s, D rs , y, Dy) 

3=0 " J 



20 



STEVE ZELDITCH 



where 

(t)nl J (- S ,D s ,y,Dy)=n^ j n 2 +n^TZ+TZ^,wzthTZ^ € C~({0,L],£r 2k ); 

z 2 2 "2 2 2 

K_ 

(ii) TZ 2 2 _ j (s,D s ,y,Dy)\ — TZ^Z°{s,y,Dy)\ = fj(I) for certain polynomials fj 
of degree j+2 on R; 

K_ 

(Hi) U 2 2k+1 (s, D s , y,D v )\ = TZ^ ,0 2k+1 (s, y, D v )\ = 0; 

Z 2 Z 2 

(c) satisfies the boundary condition: CW^ K (0)C — W^ K (0), CW^(L)C = 

" 2 

(d) We have 

[(W+ § )-^ivW+ f ]|o ~ -A^ 4 + 2N 2 TZ + F 2 NK {I) + E n k 
where F^ K (I) := J2f=o N ~ 2 Vj(-0 an d where 

(e) The error term E N+ K_(s,y,ff) = N~( k+1 ^E n k_ where I~( K+1 ^E N i£ is a bounded 
operator on L 2 (M.) for each s with a uniform bound in N . 

Proof 

After the linearization step, the scaled Laplacian has the perturbative form 

oo 

(42) TZ+ = iV 4 + N 2 TZ + E+ , E+ Q := ]T N^ m ^^ 

m=2 

where 

(43) U 2 -ni e C°°(K,£ e m ~ 4 )^ 2 + C co (R,£™- 2 )K + C cc {R,£™) 



Actually, as discussed in (2.7), a 7 fc depends only on the class of A in 5 ,2 > 2 ( fc + 2 ) (T*U, R + j) 
and by (3.4) this is the same as the class of the complete symbol a^ N of A^v at 
(7 = modulo N~ 2 ( k+2S> . We may therefore drop higher order terms to get a poly- 
nomial partial differential operator A 2fe+2 with the same wave invariants a 1 j for 
j < k. . For simplicity we do not indicate this truncation in our notation, but the 
reader is invited to think of A as a polynomial differential operator in the (y, Dy) 
variables. 

We now construct Weyl symbols Pj/2,Qj/2 so that iterated composition with 
e N 1 ^ p+l< ^h/2 w in successively remove the lower order terms in TZj^ after restric- 
tion by \ and so that the boundary condition is satisfied. As mentioned above, all 
operators will be standard Weyl pseudodifferential operators defined in a neighbor- 
hood of f2n and acting only on the transverse space M. y with coefficients in s. Hence 
we may construct them using the symbolic calculus. 

Let us first rewrite the boundary condition in terms of Weyl symbols: On the 
symbol level, we have 

(44) Ca w (s,y,D v )C = a w (s,y,~Dy) 

where a is the complex conjugate of a. Hence it is natural to split up a Weyl 
symbol into its real/imaginary parts and into its even/odd parts with respect to 
the canonical involution (y, fj) — > {y, —fj). We note that this splitting is invariantly 
defined since a is real if and only if a w (y, Dy) is self-adjoint and since the even/odd 
parts of a real symbol are its even/odd parts under conjugation by C. Since our 
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symbols are always Weyl symbols in this section, we omit the superscript w in the 
future. 

By assumption Pj/2 an d Qj/2 ar e real symbols so it remains to split them into 
their even/odd parts: 

(45) Pj/2 = P/ /a + P° /2 , Qj/2 = Q)/2 + Qj/2 ■ 

Now we observe that for general Weyl pseudodifferential operators P, Q, 

C(P + iQ)C = iQ iff P° = Q e = 0. 
Therefore, the boundary condition on Wjt is equivalent to: 

(46) P° /2 (0,y,fj)=P° /2 (L,y,fj) = 0, Q« /a (0, y, fj) = Q e j/2 (L, y, fj) = 0. 

We emphasize that there is no condition on Q°/ 2 or Pj/2- 

To begin the induction, let us construct Pi (s, y, Dy), Qx(s,y,D$) e C°°([0, L])<g> 
£f so that the boundary conditions are satisfied and so that 

(47) e-^ 1 ^*^" 1 ^* \ = {N± + N*Tl + EUlo- 
Expanding the exponential, we get to leading order the homological equation: 

(48) {[n,{p + iQ) h ] + n h }\ = Q. 

Taking the complete symbol of both sides we get the symbolic homological equation: 

(49) i{a+-I,Pi +iQ±} + 1li\ = Q. 

]_j 2 2 2 

The equation may be rewritten in the form: 

(50) d- s (P + iQ)i{s,r a (s)(y,f})) = -iKi\ . 
The solution has the form: 



(51) (P + iQ) h (5,r a (s)(y,fj)) = (P + iQ)i(0)-iJ Q Ki(u,y,f))du. 

We need to determine (P + iQ) i (0) so that the boundary conditions at s = and 
s = L are satisfied. To clarify the boundary conditions and their solvability, we 
separate out real/imaginary and even/odd parts in the equations to get: 

dgPt + p?} = | D d- s p? + f {/, Pf} = qui \ Q 

2 2 2 2 2 2 

(52) 

dsQi + f {/, Q1 } = -fRH\ \ d g Q1 + f {/, Qi } = -mil \ 

2 2 2 2 2 2 

together with the boundary conditions 

(53) Pf(0) = Pf(i) = 0, Q e 1 (0) = Ql(L) = 0. 



We now observe as in the boundaryless case (cf. Z.l], Lemma (2.22)) that 
(u, y, fj) = yo(ar ( ^ 1 (u))/ is a polynomial of degree 3 in (y, rf) in which every term 

is of odd degree in (y, rj). By (|5(]) it follows that Pi,Qi are also odd polynomials 
of degree 3. 
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To analyse the equations further we change coordinates. In the elliptic case, we 
use the complex cotangent variables z = y + ifj, z = y—ifj which satisfy {/, z m z n } = 
i(zJ^ — z-g=)z m z n — Urn — n)z m z n . We then write: 



Pf(s, 

2 


z, 


z) 




+ z m z"), 


2 


z, 


z) 




- Z m Z n ). 


Ql (s. 

2 


,z. 




= Sm,n:m+n<3 9| m „ (■*) ( zI ™^™ 


+ Z m Z n ), 


Q\ (s. 

2 


,z. 




= Em,n:mfn<3 l\mn W i 2 ™^™ 


- z m z n ). 



In the hyperbolic case, we first use real coordinates (y, rf) in which I h — \{rj 2 — y 2 )- 
We then define w = y + r],w = y — r] so that I h = ^\w\ 2 :— ^ww and so that 
{I h ,w m w n } = (m — n)y m rj a . Since the elliptic and hyperbolic cases are quite 
similar, we only carry out the details in the elliptic case and refer to |Z.2fl (Lemma 
3.1) for complete details on the hyperbolic case. We will refer to to = n terms as the 
'diagonal terms'. They only occur for even j and are terms which are polynomials 
in the action variables. 

Then ( p2| ) may be rewritten in terms of these coordinates. In the elliptic case, 
we have: 



(55) 





M§(to- 




3^! | 

2 


iip| mn ( s ) - 


Mf(m- 


n)p e i (s) = 


2 


#■£?! (s) - 


hif (m- 




2 


4-ql (s) - 


hif(m- 




-awe? | 

2 



Similarly in the hyperbolic case, although there is no factor of i in the second terms 
on the left sides. Since to ^ n, we can immediately solve for p\ ,q°i : 



v\ (s) = r L A ^P°i ( s ) ~ \o} 

(56) 

We thus eliminate the p\ , q° variables and reduce to (uncoupled) second order 
equations for the independent variable P°i mn i 9| mll ( s )- I n the elliptic case, they read: 

-&P% mn i S )- [f(™-*l mn (*) = ^y^||o+^|| 

(57) 
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In the hyperbolic case, the operator on the left side is — 4-s + [j-{ m ~ n )] 2 ■ The 
boundary conditions on P°i mn , Q± mn are ( m both elliptic and hyperbolic cases): 

(58) 

<4„ m (o) = o, q mn (L) = o. 

The boundary value problem ( |57| ) - ( p8| ) is always solvable unless is an eigenvalue 
of the operator D 2 — [-g (m — nff (elliptic case), resp. D 2 + [j^(m~ n)] 2 (hyperbolic 
case) with boundary conditions q(Q) = = q(L) The eigenfunction would have 
to have the form sin(|j(m — n)s) in the elliptic case or sinh(j^(m — n)s) in the 
hyperbolic case. Hence in the elliptic sufficient condition for solvability is 

that a/n ^ Q while in the hyperbolic case there is no obstruction if a ^ 0. 

Thus we have solved the conjugation problem to third order. The exponents 
Pi,Qi are odd polynomial differential operators of degree 3 with smooth coeffi- 
cients defined in a neighborhood of [0, L}. By construction, the Qj/2's will always 
have the same order, same order of vanishing, and same parity as the restriction 

1Zi\ . It follows that 
2 1 

N^addP + iQ)!) : N^ k ^ l (M.)(^£'" 1 -> A^ (fc+1) [^(R) ®£ e m+3 + ^(R)<g)£ £ m+1 ]. 

We now carry the process forward one more step because, as in the boundaryless 
case, the even steps require something new. In the second step, 1Z m is replaced by 

111 - e^ 1{P+lQ) hn N e N ' 1{P+lQ) h e x R). 

We expand in powers of N to get: 

00 'j 
(59) Hl~Y. N ~ A+n E j{ad{P + iQ) h yil 2 -m. 

n—o j'-f-m— n 

00 

:= A^ 4 + N- 2 K + Ar 4+ ™7£|_„. 

n—3 



An obvious induction as in [Z.l| gives that 



ad(P + iQ)iyn 2 -^ € C°°(R,£™ + i- 4 )1l 2 + C°°(R, £™ +J - 2 )1l + C°°(R, 8™ +j ). 
1 

It follows that HZ n has the same filtered structure as 
1 

We now conjugate with an exponential of the form e N ( p +' I( 2)i. As above, 
it leads to a boundary problem analogous to the previous case with j = 1 except 
that now diagonal terms with m = n do occur. In the elliptic (resp. hyperbolic) 
case, a diagonal term is a function of \z\ 2 (resp. |u>| 2 ), hence is even under the 
involution z — > z (resp. w — * w). Since the equations for the even/odd coefficients 
decouple completely, the boundary condition reduces to: 

Pf(0, \z\ 2 ) = Pf(L, \z\ 2 ) = 

(60) 

Qi(0,\z\ 2 ) = Qi(L,\z\ 2 ) = 0. 
Similarly in the hyperbolic case with w in place of z. 
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The same is true for any even j so let us consider the general elliptic case. We 
write the diagonal terms in the form: 



(61) 



p?(sM 2 )=j: m : m < j pUs)\ z \ 2m 



1 2m 



It is obviously impossible to solve the boundary problem: 
(62) 

£qf m (s) = -XR4_ S \ 

with zero boundary conditions in general. To satisfy the boundary condition we 
need to add terms /j(|z| 2 ) to the right side. Then we get: 



(63) f s pj m {s) = XRr,,. 3f,(M 2 ), £qf m (s) = -WR^\ + Ztf^z] 2 ) 
with the boundary condition above on p^ m , qj m - We can solve the equations with 



(64) 

where 
(65) 



(s) = r Au. |z| 2 ) + QfjQzftydu, 



/i(N 2 ) = -i / Ki%(u,\z\ 2 )du. 



In the case j = 1, we have then conjugated^ to N i +N 2 K+Op w (f 1 (I))+0(N- 1 ). 
We note that Op w (fi(I)) is a function of /, so we have conjugated to normal form 
to fourth order. 

We then proceed inductively to define polynomial symbols (P + iQ) ± (s, y , fj) , 
polynomials fj (I) and unitary iV-pseudodifferential operators W N i := exp(N~i (P+ 
iQ)i) such that: 

d) ni 

<2.K 



N 

a) n 



N 



w-%w-^...w-Itz nWj 



N 4 + N 2 K + J2 n = 3 N 4 - n K 2 2 _ 



W N K-i_W N K\ 



a i) 



for K > n-2,1lZ „ =K 



, ,v 2 _„ - ,V2_» = [K, (P + iQ)n_i] + n'_ 



iv) n^&D-^D^o = K™_l°(s,y,Dy)\ = /,(/)| 

v) 2K+1 (s,D s ,y,Dy)\ ^TZ^'° 2K+1 (s,y,Dy)\ = 0; 

vi) K 2 2 _™ e C°°(IR,£ £ m " 4 )^ 2 + C°°(R,£™- 2 )K + C°° 

vii) {P + iQ)j< e C 00 (R,£7 +1 ); 



cm \ . 

1 C £ )1 



(viii) (P° + iQ e ) k(0) = (P° + iQ e )K(L) =0; 

(ix) E+ K is divisible by jV 4 -( 2ir+1 ). * 

" 2 

The details about the degrees and parities of the polynomials are similar to the 
boundaryless case of [Z.l| and we therefore omit the details. An important point 
which wc must address here is that the real parts Pj/2 do not contribute to the 
principal part of the exponent (P + iQ)j/2- More precisely, the homogeneous part 
of Pj/2 of leading order 2j + 2 equals zero. This happens essentially because the 
Weyl symbols of A, An and 1Z N are real, all three operators being self-adjoint with 
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respect to Lebesgue measure. Conjugation by unitary operators e lN 3 Qi/ 2 preserves 
this reality. At first sight it appears that the 72.|_ m operators continue to have real 
symbols and therefore that Pj/2 is zero. However, the | operation kills part of 
the operator and in particular it kills the self-adjointness and reality of symbols. 

i 

The 7£|_m generally have complex symbols and one does needs to conjugate with 

e N- j (P j/2 +tQ j/2 )_ 

What is true is that the leading order terms of order m the Weyl symbol of 

k 

7?.|_m in (y,r]) is real. As is visible from (vi) aboves, such terms have no factor of 
1Z in front while any term with a factor of 1Z has a lower degree in (y, 77). Hence 
taking \ a commutes with setting an = in terms of top degree. But setting an = 
in a real symbol produces another real symbol, so the terms of top degree in (y, rj) 
must be real. 

Let us consider the statement (ix) about the error terms. In each conjugation we 
expand the exponential up to order N~ K and absorb the remainder in the higher 
powers of A -1 . 

To analyse the error, we write the partial Taylor expansion with remainder of 
the exponential function as: 

e lx = e M (ix) + r M (ix), e M (ix) = 1 + ix H h 

r M (ix) = (ix) M + l b M {ix), b M {ix) = £■ ■ ■ J* t%t™z\ ■ • • t° e t -' t '--^ t °™ dt ■ ■ ■ dt M - 

Clearly bu is a bounded function on the axis of real x. Now plug in ix — A~- J ad(P+ 
iQ)j/2)', although it is complex, its principal part is real. We get that 

e N-^d(P+iQ) j/2 = j + N -i ad rp + iQ) + . . . + (N-^d(P+iQ) j/2 )r 

(66) 

+[N-Jad(P + iQ) j/2 ] M+1 b M (N-^d{P + iQ) j/2 ). 

At the 2ATth stage we are conjugating with 2K factors e lN 3 ® j/2 ,j — 1, . . . , K. 
To get a normal form up to order N /L ~^ 2K+1 ^ at the Kth stage it suffices to choose 
each Mj such that j(Mj + 1) > K. The eja terms give a polynomial differential 
operator satisfying (vi) and the bu terms give the error E~t K . By construction, it 

is N 4 ~( 2k+lS> times a conjugate of a polynomial differential operator of degree K + 2 
in (37,17). □ 

4.4.1. Remarks, (a) The operators WjT with which we began are now determined 
by ©, i.e. by W+ = T Nk W+T^. We have: 

(67) W+ = n°° =1 e N ~ jTN >° {P+iQ)i/2T ^. 

As will be clear below the Wk's belong to the class of homogeneous Fourier integral 
operators, unlike the W^s, which are exponentials of isotropic pseudodifferential 
operators. 

(b) In the next section we will glue together the Wv"Vs into a Fourier integral 
operator. The vanishing of the principal part of Pj/2 then implies that it does not 
contribute to the principal symbol of (P + iQ)j/2, but only to the amplitude. It 
will also imply that &M(A~ J ad(P + iQ)j/ 2 ) is a bounded FIO. 
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(c) Finally, let us note that CW+e lN * L C = W+e m2L if and only if e- lN * L = e lN * L , 
i.e. if and only if N 2 — jk for some k. Henceforth we put N k = y/j;k and only 
consider these special values of N. For notationally simplicity we write Wjt for 
W^f . We also define F(k,I) as the formal asymptotic series whose -ftTth partial 
sum is F NkK (I). 

5. Conjugation of the Dirichlet Wave group 

So far, we have conjugated the semiclassical Laplacian on the open space in a 
microlocal neighborhood of 7 to a semiclassical normal form modulo a small semi- 
classical remainder. We now glue together the component intertwining operators 
W k = TN k W£T^ k into a homogeneous Fourier integral intertwining operator of 
the Dirichlet wave group to its normal form. We emphasize that the intertwining 
operators are defined in the open space. 

Conjugation by the scaling operator T and its local form Tjv fc plays an important 
role since it converts isotropic pseudodifferential operators into standard homoge- 
neous pseudodifferential operators. This will be discussed further below. Since it 
becomes heavy notationally to distinguish the action operators I from their con- 
jugates we abuse notation somewhat by denoting both I and TIT^ 1 by the same 
symbol I. 

By the Dirichlet wave kernel we mean the fundamental solution E(t,x,y) of 
the mixed wave equation (^) with Dirichlet boundary conditions, i.e. the kernel 
of cost\/ An). The 'free' Laplacian A on the open space is only well-defined in a 
neighborhood of 7, so to be precise we need to cut it off to U e . We will not indicate 
this cutoff in the notation because we will explicitly microlocalize it later on. 

5.1. Local components on S\ L x HL By 'local component' we mean an operator 
which only acts on a specific Fourier coefficient. 

Definition 5.1. Define the operators: W k : L 2 (S\ L xR)^ L 2 (S\ L x R) 

W k = T e {W+U +k - W-Il- k ). 

Here, W k ~ = CW+U +k C. 
In 

Proposition 5.2. We have: 
ft) 

r e AVF±n ±fc = Te^FdDl, J) 2 II ±fc + W^E^U^}. 

Here, \D\ = y/-D*D on S\ L . 
(ii) 

r e Aw k u = T e {w k F(\D\,i) 2 n + w k E+ K n }. 

(Hi) The Schwarz kernel ofW k Tl vanishes on dQ . 



Proof 

(i) This is essentially just a restatement of Lemma (4.6). We can prove it by testing 
both sides against functions of the form e lks f(N k y). We do the case k > 0. On the 
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set r £ = 1 we have 

^ ft AT. „,\ — ^i^^-i 2 ^ A^i^T 1 .. M, 

k N. 



AW+e^f(N k y) = e^e^Hr Ae^T Nk W+T^f(N k y) 



= e^T Nk A Nk W+f(y) = e i2L i~T Nk W k + (F(k, I) 2 + E k K )f(y) 
= W+{F{kJ)) 2 e^f{N k y) + W+ E* f (N k y) 

= {W+(F(\D\J)) 2 + W+E* )}U ±k e^f(N k y). 

Between lines two and three, the action operator I changed from its isotropic form 
to its homogeneous form. The case of k < follows by taking complex conjugates. 

(ii) This follows immediately from (i) by taking the difference of the + and — cases. 

(iii) The Schwarz kernel of WH is given by '^2 lk W k Ii ok , so it suffices to show 
that W sin(^p) / (N k y) = on ,s = 0, L. This follows from Lemma (p~6|). since 

(W k sm(^) f(N k y) = on s = 0, L. Also, we obviously have sin(^-)f(N k y') = 
on s = 0, L. □ 

5.2. The homogeneous intertwiner. We glue the W k together as follows: 

Definition 5.3. Define the operators W ± ,W : L 2 (S\ L x Rg) -> L 2 {S\ L x E s ) by 

(a) ^- = Er=i^n- 
(iii) w = w+ -w~. 

We now analyse the Fourier integral nature of W, W . A key point is that the 
operators T^ k glue together to form a 'global' transverse scaling operator 

T : L 2 ([0,L} s x K s ) -> L 2 ([0,L} s x %) 

(68) 

T S m(^)f(y) :=N kS in(^)f(N k y). 

Conjugation with T transforms the isotropic Weyl calculus into the usual homoge- 
neous psuedodifferential calculus. For instance 

(69) 

T(D 2 + f)T* = \D s \- l D 2 + f\D s \, T(D 2 - f)T* = \D s \- l D 2 - f\D s \. 

As is verified in |g[ and elsewhere, T is an oscillatory integral operator associated 
to the canonical transformation 

(70) tp(s,a,y,n) = (s + ^fay.^fa 1 n). 



The following proposition is analogous way to (|Z.l], Proposition (3.4)). Below, 
the notation A ~ B in V means that A — B is smoothing in V. In the following 
proposition we assume the order K of the Birkhoff normal form is infinity. 

Proposition 5.4. There exist conic neighborhoods V ± , V ± of := 7 C H {±er > 
0} and canonical transformations \ ± : V ± — > V ± such that: 

(i) y± - Id on M 7 ± := {a, a,y,fj)}&V ± :y = f] = 0}; 

(ii) W ± e I n (Sl L x R x S\ L x R, grx ± ) where grx^ C V ± x V*' is the graph of 

x ± ; 
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(Hi) W ± is elliptic in V ± x V ± . 
Proof 

To analyse the sum over k, we enlarge the 'open model space' S\ L x i to the 
'space-time' S\ L x R x S\ L and define the operator 

(71) 

e it\D s \ . L 2( S i L x M ) ^ L 2( S i L x | x sl L ), e u \ D ^e^f(y) = f{y). 

The range of e**' 135 ' is contained in the kernel Ti of the 'wave operator' |-D S | 2 — \D t \ 2 
on S\ L x M x S\ L . We also denote by P the orthogonal projection to 7i. We then 
introduce the further operators: 

Definition 5.5. Define W^(s, y, D v ) and W\ Dt \ on L?(S\ L x R x S\ L ) 

(a) For k > put 

W+ tl (s,y,Dy)Il + f(s,N k y)e i2 ¥ = e^W+ f(s, N k y); 

For k < put it equal to zero. 

(b) W lDtl (s, y, D,) = CW+ t] (s, y, D y )C; 

(c) W\ Dt \(a, y, D v ) = W+ tl (s, y, £> s )n+ - W^{s, y, D y )IL_ . 

These definitions suggest the introduction of a scaling operator T adapted to 
the ^-variable. We therefore define fe ltm f(s,y) = N m e ztm /(s, N m y). We are only 
interested in its action on the invariant subspace Ti. where the frequency in s and 
t are the same. As with T, T is an oscillatory integral operator with underlying 
canonical transformation 

i>(s,o-,t,T,y,r}) = (s,a,t + — ,t, V^y, sJt n). 

We then have: 

W+ Dtl (s,y,D y ) = U^ Dt ^ 3/2< - p+ ^^'^ s -y' D y\ with P = TPf-\ Q = TQT' 1 . 

The operators P, Q are (usual) homogeneous pseudodifferential operators, as we 
will argue below. 

The following identity is the key to the Fourier integral properties of W ± ,W: 
Let j : S\ L x R — > S\ L x S\ L x M be the inclusion j(s,y) = (s,0,y). Then: 

W ± {s, y, D v ) = j*W± tl (s, y, Dy)e*\ D '\ 

(72) 

W = j*W\ Dtl e^ D ^ 

To prove it, we apply both sides to functions of the form f(N k y). In the plus 
case, we have (for k > 0), 

j'W+ t \(8,y,Dy)e*\ D '\e iS Pf(N k y) 
= j*e^W k + (s,y,Dy)e^ f(N k y) 

= W+e^f(N k y) 
and similarly in the minus case. 
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We now complete the proof of (i)-(iv). For simplicity we only consider the 
+ case, the other being essentially the same. By definition, is a product of 

factorsof the form e \ Dt \ ( p i/* +iQ i/*) w j t h P j/2 ,Q j/2 scaled polynomial differential 
operators of degree j + 2 in the variables \D t \^y and \D t \~?D y . As noted in the 
proof of Lemma (|4.6| ), the terms in Pj/2 have the parity of j and the leading 
order term in (y,D y ) vanishes. Hence, \D t \~^ 2 Pj/2 is actually a pseudodifferential 
operator of at most zero order and its does not affect the status of WX )t , as a Fourier 
integral operator. Since Qj/% is a polynomial order j + 2 in (\D t |-Dt| — ^D y ), it 
follows that | A|" J ' /2 (P + / 2 is a first order (homogeneous) pseudodifferential 
operator of real principal type for each j. Its exponential e l ' Dt ' J/ is 
therefore a Fourier integral operator. Hence the product of any finite number 
of factors is a Fourier integral operator. As discussed above, we only need to 
deal with finitely many factors, but we note that the the full infinite product can 
be regularized as a Fourier integral operator. Indeed, one may apply the Baker- 
Campbell-Haussdorf formula to any finite number of factors to rewrite | in the 

form W7j)A = e^j =1 ' £>t ' 3 ' f p i/ 2 +'2i/2)_ jf onc interprets the sum as given by the 
Borel summation method then the exponent is well-defined as a pseudodifferential 
operator of order one and of real principal type. The principal symbol of the series 
J2j \Dt\~^ 2 {V + iQ)j/2 is then well-defined as a an element of S 1 ' 2 of the form 

H(s,y,V: T ) : = Ej t~ j/2 9j/2(s, VnA in (y,v)i witn Qj/2 tne leading order 

homogeneous part of the Weyl symbol of Qj/2- The symbol of Vj/2 is of order at 
most zero as a symbol and hence does not contribute to H. 

The operators e 1 *'- 03 ' and j* are clearly Fourier integral operators, with associ- 
ated canonical relations 

C={((s,a,y,rj);(s + t,a, t, a, y, r,))} C T*{S l 2L x R x S\ L x S\ L x E) 

(73) 

I> ={((s,a,y,ri);(s,a,0,a,y,r))} CT«(4 xRx^ xRxS$ L ). 

The canonical relation underlying j*W^ D ^e' lt ^ D ^ is therefore given by the com- 
posite relation 

A = r> o gr (j) 1 oC 

where <j) u — (expwS^) with 3jj the Hamiton vector field of H and expu2# its flow. 
The first two factors compose to the relation {((s, a, y, rj)\ (^{s -I- 1, cr, t, a, y, i]))}. It 
is easy to see that ^(s + t, cr, t, a, y, rf) has the form (s + t, cr(l), t(l), cr, y(l), T](l)) 
where x(l) stands for the value at u = 1 of the of x-coordinate of the orbit of 
<fi u through the initial point (s + t, cr, t, cr, ?/, rj). To contribute to A one must have 
= t(l) = t + J Q Q^r(4> u (s, cr, t, cr, y, r])du. We observe that H vanishes to order at 
least two along R + j = {y = rj = 0} and hence that <j> u acts as the identity on 
this set. In particular, t(l) =t = Q there. Futhermore, for (y,rj) sufficiently small, 
the derivative of i(l) with respect to t is non-zero and hence there exists a unique 
solution t(s, cr, y, rj) of the equation relating i(l) and t. It follows that (at least) 
in a sufficiently small cone around {y = rj = 0}, A is the graph of the canonical 
transformation x + (s, cr, y, ?/) = (p (s + t(s,<T,y,rj),<T,t(s,<T,y,r)),(r,y,r]) (where the 
(t, r)-coordinates are omitted on the left side). 

The — component has a similar description. Therefore, W is a microlocal Fourier 
integral operator associated to the graph of a canonical transformation defined in 
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a small cone around M + 7 . Each of the three operators in its composition has a 
canonical principal symbol and it follows in a standard way that o~(W) is a graph 
1/2-density. Hence W is microlocally elliptic, concluding the proof. □ 
Below we will need a somewhat smaller open cone with the following property: 

Proposition 5.6. Let int(-f ) = -f Q n T*(intQ ) where intQ Q denotes the interior 
ofil . Then there exists an open conic neighborhood V ofint(j ) with the property 
that x(Vo) H es i n the interior o/T*(f2 Q ). 

Proof The point we must address is that \ could carry a covector at (s, y) G int(fl ) 
to a covector on <9f2 Q or in the exterior. To determine whether this happens 
we must study s' — -§^H. We observe that this component is homogeneous of 
degree on T*fl ) and that it vanishes to order three along y = rj = 0. It 
follows that there is a neighborhood V a of -f a of the form max{|s|,|L — s\} < 

Cmm{\y\ 3 ,\y\ 2 ^,\y\ 1 -^-, 1 -^-} with the property that <j>(s,y,o-,r)) e T*(mift ) if 

(s,y,v,v) e V a . □ 

5.3. The error term. We now make a similar analysis of the error term. We first 
define operators E ±K on fl±(S| L x M) and E K on L 2 {S\ L x E) by 

E+e ik ^f(N k y) = W+E+ kK e ik ? f(N k y), (k > 0) 

E- K e ik ^f(N k y) = CW+E+ kK Ce ik ^f(N k y), (k < 0) 

E K = E+U+ - E~U_ 
To analyse the kernels of E^ and Ek we introduce the operator on S\ L xlx S\ L : 
(74) E+ Dt{K = W+ D ] l AW+ Dtl - Fl(\D\J). 

Recall here that F 2 K {\D\J) = Y.)=i fj(\D t \J)- 

Proposition 5.7. As an operator on L 2 (Sl L xR), we have: E~^ = j*Wy D ^E^ D ^ K e 1 ^ 
Proof 

As above, we test both sides on functions of the form e lk ^ f(N k y) for k > 0: 

= j*W+E+ kK e i ^e ik ^f(N k y) 

= W+E+ kK e ik ^f(N k y) = E+e ik ? f(N k y). 
Similarly for k < 0. 

□ 

In the following we denote by V the cone in T*(S\ L x E x S\ L ) defined by 
(s,a,y,7],t,T) e V iff (s,a,y,rf) G V,C < \r\/a < l/C for some C < 1. We also 
denote by R + 7 x T* S\ L the symplectic symplectic subcone of V which is defined 
by y = V = 0. 

Proposition 5.8. E+ Dt]K e OpS 2 ' K (V, E+ 7 x T*S$ L ). 
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Proof By (]74|), by Proposition (5.6) and by Egorov's theorem, we see that E+ K 
is a homogeneous pseudodiffcrcntial operator of order 2. To show that its symbol 



lies in S 2 ' K it suffices by proposition (15) to show that the Taylor expansion of its 
scaled symbol is divisible by jV - C K '+ 1 ). The symbol of £t , K equals E~^ K with |r| 
substituted for jV 2 and we will denote it by E^ K (s, y, rf). Here we use that E^ K 

is independent of a. Rescaling in V replaces |r| by TV 2 . Moreover, the transversal 
scaling (jfU, Nrj) is canceled by the operator T. Hence the scaled symbol is precisely 
the original semiclassical symbol ENK(s,y,fj)- By lemma ( |4.6D it is divisible by 

N -(K+1). ■ q 

From the previous two propositions we have: 

Corollary 5.9. E/<- is a sum of terms of the form ARB where A,B are bounded 
Fourier integral operators and where R £ OpS 2 ' K (V, IR + 7). 

5.4. Conclusion. The following proposition sums up the discussion. It gives the 



homogeneous analogue of Proposition (5J2) and the proof is essentially the same. 



We use the notation A ~ B in V to mean that A — B is a Fourier integral operator 
of order — oo in V . Also we emphasize that the action variables are the scaled ones, 
i.e. functions of |-D|2y and \D\~^D y . 

Proposition 5.10. We have: 

(i) AW ~ W[F(\D\J) 2 +F, K )} inV. 
(it) Ifue ffo(^o) then Wu = on dn o . 

(Hi) The canonical transformation \ conjugates the Hamiltonian on V to 

the Hamiltonian F(a, I) on Vj 1 modulo an error which vanishes to infinite order 
on R+j . 

Proof 



(i) Since W = Y,k( W k U +k ~ ^ fe ~ n -fe) it follows as in Proposition (|J) that 



T e AWe ik ^f(N k y) = T e W[F(\D\J) 2 + E K ]e^ f(N k y). 

(ii) This follows by taking the principal symbol of the equation in (i). In the 
principal symbol, the remainder is homogeneous of order 1 and vanishes to order 
K at 7 . 

(iii) It suffices to show that W sin k^f(Nky) = on dft . But this follows from 
the fact that W k smk^f(N k y) = on dn o for all k. 

□ 

5.5. Normal form and microlocal parametrix. We now use W to conjugate 
the (odd part of the) wave group of the normal form to a kind of parametrix for 
the mixed problem on Q Q . First, let us define the odd Dirichlct normal form wave 
group: 

Definition 5.11. By F (t, x, y) we denote the odd part of the fundamental solution 
of the wave equation: 

- F(\D\J) 2 F (t) onRx {S\ L x E) x (S^ L x R) 



at 2 



f o (0) = n G (o) = o 



F o {t,;z)=0 = F o (t,z r ) on dfl, 
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To be correct, F is only denned in a microlocal neighborhood of j . We do not 
indicate this in the notation since we will microlocalize it later on. In the elliptic 
case, for instance, F a has the eigenfunction expansion 

(75) F (t, a, y; s', y') = ^ cost(F(fc, q + y)cj> kq (s', y') 

k,q 

where 4>k q {s,y) = sin k^D q (Nky). By odd part we refer to the projection under 

n . 

We now wish to conjugate the odd normal form wave group to a microlocal 
paramctrix for the Dirichlet wave group E(t). 

Proposition 5.12. There exists a conic neighborhood V' C V of j and a zeroth 
order pseudodifferential operator G such that: 

WW*G ~ W*GW = I + R with WF[R) n V = 

WF 2 W*G = A + A 1 E K A 2 + Ri, WF(Ri) n V = 
where A\,A 2 are zeroth order Fourier integral operators. 



Proof: By proposition (5.4), W is a zeroth order Fourier integral operator micro- 
supported iiifx x(V) an d associated to the graph of the canonical transformation 
X- Hence WW* is a zeroth order pseudodifferential operator microsupported in V 
with symbol identically equal to one in a smaller cone V C V. Hence there exists a 
positive zeroth order self-adjoint pseudodifferential operator G microsupported in 
V such that WW*GI + R with WF{R) n V = 0. 
Let us put W-i — WG. Then by proposition (5.10), 

WF 2 W-i = WF 2 W*G 
(76) - AWW*G + WE K W*G in V 

- A + R in V 

where R = W~Ei K W*G. □ 
We now introduce a kind of microlocal parametrix for the wave kernel. First we 
define a microlocal cutoff ip to the conic neighborhood V of M.^ a , with tp = 1 in a 
conic neighborhood of intW + j. 

Definition 5.13. The parametrix is defined by: £(t) = WF T(t)ipv^n W-i. 

The extent to which £(t) is indeed a microlocal parametrix near 7 Q is given in 
the following lemma: 

Lemma 5.1 4. F or any A 6 5 , *(J7) with microsupport WF(A) C V (defined in 
proposition ( \5.dj) , E(t)A and £{t)A are Fourier integral operators in the class 
I°{{Sl L x M) x [S\ L x R), grx) + / ((^ L x M) x {S\ L x M), gr\ or). We have 

-E(t)A~£(t)A + R(t)A 

where R(t) is a finite sum (or smooth compactly supported integrals) of terms 
B\Ri{t)B2 where Bi,B 2 are bounded Fourier integral operators and R\ € OpS 2,K . 
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Proof: E(t) is uniquely characterized as the microlocal solution of the Cauchy 
problem (^) on fl a with initial condition equal to the identity operator on L 2 (Vt ). 
Hence it suffices to show that £ is a also a microlocal solution modulo errors in the 
stated class. 

We first verify that £t) is a microlocal solution of the wave equation in V'. But 
(77) 

^£{t) = WF(\D\,i) 2 F Tp v l Qo W-i 

= WF(\D\J) 2 W-iWF rP v ln W-i+WF(\D\J) 2 (I-W-iW)F i) V l Qo W- 1 
= A£(t)+E K £(t)+R 1 (t), 

where Ri (t) = WF(\D\J) 2 {I-W- 1 W)F ^ v l Qo W- 1 +R 2 £(t) with WF(R 2 )nV = 
0. By shrinking V if necessary, the first term also has order — oo in V. Hence £{t) 
is a microlocal solution of the forced wave equation with forcing term of the form 
E K £ t {t)+R 1 {t). 

Regarding the boundary condition, we observe that WH (x,y) = if x € dft 
and so £ (x, y) = if x £ dQ. 

The initial condition is more difficult. We must show that WH ipv^n W-iA ~ 
^^4 for a pseudodifferential operator A microsupported in V To this end we first 
recall that n o = J (J — r) where r is the reflection (s, y) — > (— s, y) (with s mod 2L.) 
Hence WH ipvln W-iA = \W Vv l n „ - ±W r lo^WLiA It suffices 

to show that WiJ> v ln W-iA ~ A and that W r ijj v ln o W-iA ~ 0. 

The operator Wipv^n W-\ is a pseudodifferential operator with a singular sym- 
bol. Temporarily ignoring the singularity along <9f2 we may apply Egorov's the- 
orem formally to find that the complete symbol is identically equal to one in V . 
To see this, it is convenient (although not necessary) to use that W = e fo r 
some first order pseudodifferential P of real principal type (see proposition ( |5.4| )) 
We then put W(u) = e luP and consider the conjugation W '(u)ipyln W-i(u) . In 
a well-known way (cf. §7.8), the complete symbol expansion of this opera- 

tor is obtained recursively by solving transport equations along the orbits of the 
Hamilton flow (j> u of dp = H. The principal symbol equals <p u * (i/v1q o ). Simi- 
larly the complete symbol at (/^(x,^) of Wipv^n W-i depends only on the germ 
of the complete symbol of tpv^n at (^(x, £) and on the germ of (j) 1 at (x, £). Since 
1 (V O ) C V C\ T*(intSl ), the germ of the complete symbol of tpvlQ is identically 
one at (a;, £) G V and hence the germ of the complete symbol of Wipvlsi W-i 
equals that of WW-i at <j>(x,£), hence equals one identically. Since we have com- 
posed with A satisfying WF'(A) C V this calculation is valid on microsupport of 
A and thus Wip v hi W- 1 A ~ A. 

The statement W r i/)yln o W-iA ~ follows from the fact that the underly- 
ing canonical transformation of this operator is (0 1 ) _1 r*0 1 , and the graph of this 
canonical transformation is disjoint from the T*{intVL ) n V x T*(inffi ) H V Q . This 
follows since the graph of r* does not intersect T*(intVt ) x T*(intVL ). 

Thus £(£) is a microlocal solution of the mixed Cauchy problem modulo errors in 
OpS 2,K . Since the Cauchy problem is well-posed, we have by Duhamel's principle 
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that 

(78) (S(t)-E(t)~ [ G {t-u)R{u)E{u)du 

Jo 

where R(u) is an error in the stated class and where where G is the kernel of 
smt(y f SH' ) j e Q £ ^ e mixed problem 

{df - A)G Q = on ft x n o 

(79) G (t,s,y,s',y')|t = o = &G (t, a, y, s', y')|t=o = Id 

G (t, s, y, s',y') = for s = 0, L. 
Since we are working microlocally near 7, we may replace G Q modulo smooth- 



ing operators by a Fourier integral parametrix (cf | GM | |PS[). Thus J Q G (t — 
u)R{u)E e Q) (u, s, y, s', y')du is a sum of terms of the form AR\B with R\ as above. 

□ 

6. Wave invariants and normal form 

We now prove the crucial lemma that the wave invariants of the normal form 
at 7™, i.e. the singularities of Tr cos t(F(\D\, I) at s = 2mL, agree with the wave 
invariants of y/~K at 7 m , i.e. the singularity of Tre 11 ^ at t = 2mL. Here, the 
\D\ in /) refers to the Dirichlet D on [0, L) and the I as above refers to the 

homogeneous action operator. 

Lemma 6.1. Let 7 denote a non-degenerate bouncing ball orbit of Q and let j a 
denote the corresponding bouncing ball orbit 7 Q of the normal form glow of 
Then for all m, we have a 7 ™ fe (v / A) = a 7 ™ fc (F(|_D|, I). 

Proof 

To prevent confusion between traces on [0,L] x R and traces on S 1 x R we will 
exclusively use the notation 'res' for the latter space and will explicitly put in the 
cutoff ln a (the characteristic function of f2 Q ) to indicate traces on the former. We 
have: 



a 7 fc(\/Aa) = a lok (y/ A n J 

(80) 

= res {D k t E{t)\ t = 2mL )\n . 

The residue is given by the integral of a density over the fixed point set of the 
underlying canonical relation. The density is calculated using the the method 
of stationary phase on a manifold with boundary (see Theorem ( |2.2| ) or ]GM | ) . 
essentially as in the boundaryless case. Thus the residue depends only on the 
kernel of E(t) microlocally in the cone V of proposition ( |5.6| ). Hence by lemma 
( [3.14 ) we have 

(81) 

res (D*E{t)\ t=2mL )l Qo = Ires {D\E (i)| t=2m i)ln + 2res {D k t R(t)\ t=2mL )ln 



AK A 



with R(t) G AiOpS • A 2 ). By Proposition (^6) we can drop the second term. 
Then by the definition of £{t) we get that (^l|) equals 



(82) 



2resln W{D k t F {t)\ t=2mL )^ v l no W^. 
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To evalute this residue let us approximate 1q by a smooth cutoff t c (s), equal 
to one on f2 Q and supported for s € (— e, tt + e). We have: 



a 7mfc (VA) = 2res l Qo W (D* F (t)| t=2mi )Vvln W-i 

(83) = 21im e ^ re S r e (s)W(i? t fc i? (t)|t=2 ro z^ 1 n W-i 

= 21im e ^ res ln W_iT e W(I}*F (t)| t=2mi )Vv 
where in the last line we used the tracial property of res. We note that 

W-lT e W(D*F (t)\ t =2mL)i>V 



X 



is a (standard) Fourier integral operator on the boundaryless manifold S\ L 
As mentioned above, the residue of 1q o times this operator equals the integral 
of a density over the fixed point set of the underlying canonical relation. Since 
W^\T t W is a pseudodifferential operator, the canonical relation is simply that of 
(Dj F (t)\t=2 m L) ■ Due to the factor of ln the residue density may be calculated 
as if T e = 1. But then the factors of W and W-\ cancel and we are left with 



(84) a rt (VA) = 2res l 0o (D*F a (t)\ t=2mL )^ v . 
To complete the proof we must show that this residue equals 

(85) res ln (D k t F (t)\ t=2mL )^ v . 

The difference is just that the odd kernel equals fl = ^{1 — r) composed with 
the full kernel. The two terms of n o give rise to two components to the canonical 
relation of F a . They consist of the graph of the Hamilton flow of the normal form , 
which we denote by G^™ L , and the graph of G^™ L o r. The latter map has no fixed 
points in T* (inttt ). Hence the residue of the second term equals zero. In the first 
term the factors of 2 cancel we give the result claimed in (jB5|) . □ 



Corollary 6.2. The quantum Birkhoff normal coefficients, i.e. the coefficients of 
the polynomials Pk{I) ire spectral invariants of A. 

The proof of the corollary from the lemma is identical to that in [Q] (see also 
|Z.2j| ), so we omit the proof. 

6.1. Conclusion of Proof of Theorem. We now complete the proof of the main 
result. Since the quantum Birkhoff normal coefficients of \f~K at 7 are spectral 
invariants, it follows afortiori (by taking the interior symbol) that the Birkhoff 
normal form for the metric Hamiltonian H = of the metric g$ a at the bouncing 
ball orbit is a spectral invariant. Thus we may write: 

(86) ff = H + ^ + ^) + ... + ^) + ... 

L \a\ \a\ k 

where pk is the homogeneous polynomial part of degree k + 1 in the quantum 
normal form. We now observe that this normal form of the Hamiltonian induces 
the Birkhoff normal form of the Poincare map V-y . 

Indeed, we first note that the coordinate s is dual to a and hence is an angle 
variable. Hence the Hamilton flow in action-angle variables takes the form: 

(87) <f> t (s,<l),\a\,I) = (s + tu>\ tT \,(l> + tcj J ,\a\,r), = d\ a \H,u>i = diH. 
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The billiard map from the bottom component of the boundary to the top component 
is given by: 

(88) /M<A>-0 = (<£ + *i(& W\J)vi,i) 

where t\(4>, |er|, I) is the time until the trajectory through the initial vector defined 
by (0, <j>, \a\, I) hits the upper part of the boundary. We obviously have: 

(89) h{4>,\a\,I)=L/u wv 
Similarly, the billiard map for the return trip is given by 

(90) / 8a(0,/) = (^ + t(^|<7|,/)w / ,J). 

It follows that the Poincare map has the Birhoff normal form 

(91) V 1 {4>,I) = {4> + t l ui + huiio0 ll I). 

Thus, the Birkhoff normal form of P 7 is a spectral invariant. We now want to 



apply the argument of Colin de Verdiere ([CV], §4, THEOREME) to conclude that 



the domain is determined by its spectrum. Although Colin de Verdiere makes the 
assumption that 7 is a non-degenerate elliptic orbit, his argument is equally valid for 
hyperbolic bouncing ball orbits. Let us briefly explain the necessary modifications. 

Assuming that the non-degenerate bouncing ball orbit is the vertical axis, we 
write the graph of f2 over the horizontal axis as y = f(x) = 1 +ciqx 2 + ■ ■ ■ + a n x 2n ~ 2 . 
Only even terms appear due to the left/right symmetry assumption. By definition, 
the Birkhoff normal form of the (non-linear) Poincare map P 7 is an expression for 
this map in local action-angle variables (/, 8) on the transversal. In the elliptic 
case we need to assume (as above) that a /it ^ Q to ensure that the normal form 



exists; no assumption is needed in the hyperbolic case. Following [CV| we write the 
normal form as 

T{1, 6) = (I + 0(I°°),6 + b a + b 1 I + --- + 0(1°°)). 

The key assertion is that there is an 'upper-triangular' bijection between the Taylor 
coefficients cij and the Birkhoff normal form coefficients bk- This means (i) that 
b n = B + Ca n where B, C depend only on ao, . . . , a„_i and (ii) that C^O. Granted 
(i) - (ii), {ao, . . . , &„} — > {bo, ■ ■ ■ , b n } can be inverted and the Taylor coefficients are 
determined by the normal form coefficients. 

The proof of (i)-(ii) is almost exactly the same in the non-degenerate elliptic 
and hyperbolic cases. To convince the reader of this, we briefly recall the argument 
in |CVfl and extend it to the hyperbolic case. The starting point is that one may 
write down a generating function <f> for T in terms of /: 

graph T = {(*,-;*!,-— : - = 0} 

with ()>(x, xi,s) = [(x — s) 2 + (/(s) 2 ]^ + [(xi — s) 2 + (f(s) 2 }^. Here, x, xi, s denote 
points on the horizontal axis. In either the elliptic or hyperbolic case, one may 
expand <j) = 2 + (j>Q + ■ ■ ■ + (f) n + ■ ■ ■ where (f>o (resp. (f> n ) is homogeneous of degree 2 
(resp. degree n + 2) in (x,xi,s). One easily finds that 4>q(x,xi, s) = \{x 2 + x 2 ) — 
s(x + xi) + (2ao + l)s 2 ) and that 4> n — tp n + 2a n s 2n+2 where i\) n depends only on 
ao, . . . , a„_i. The linear Poincare map P 7 has generating function 4>o and one finds 
that 

p = ( A-l -A 
7 2-4 A-l 
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where A — 2(2ao 
a > or a < — 



1). In the elliptic case ao 6 (- 
We assume A^O, i.e. a ^ 



Using the generating function one shows (|CV|, Lemma 1) that 



0) while in the hyperbolic case 



(92) T 



_ rpO 



C(x-t) 



2n+l, 



1 

-1 



+ O(\(x 1 0\ 2n+2 ) 



where T° is the (2n + 1) jet of T computed with the assumption that a n = 
and where C ^ 0. The proof is a formal manipulation with Taylor series and only 
in the evaluation of C does it matter whether the linear part of T is elliptic or 
hyperbolic. The constant C has the form C\A 2n+2 where C\ is universal. Hence it 
is non- vanishing as long as A ^ 0. 

From the construction of the Birkhoff invariants (cf. SM|) and from rtol) it 
follows that b n = B + Ca n . It remains to show that C ^ 0. To prove this, we 
recall that two germs of area-preserving transformations of (]R 2 , 0) for which the 
eigenvalues of the linear parts are not roots of unity are symplectically equivalent 



if and only if they have identical normal forms ([3M], p. 162). To bring this down 
to a finite dimensional statement, we define Gin+\ as the group of (2n+l)-jets of 
area-preserving transformations of (R 2 ,0), and I^n+i C Gan+i as the subgroup of 
elements of the form Id + 0(\{x,£)\ 2n+1 ). Also let O t = {T' € G 2n +i :T' — T = 
0(\(x,£,)\ 2n+1 )}. Assuming the eigenvalues of the linear part of T' e Ot are not 
roots of unity, the orbit r2„+i -T' C Ot consists of elements T" € Ot with the same 
Birkhoff normal form invariants up to and including b n . Put T' = T(ao, . . . , 0). If 
C = 0, then T(ao, ai, . . . , a n ) would on the orbit T2 n +i ■ T(ag, . . . , 0) for all a n , 
i.e. the curve a n — » T(oo, fflij • ■ ■ would be tangent to this orbit. However, it 



is shown in |CV| that it is transversal to the orbit. In both elliptic and hyperbolic 
cases, we can choose coordinates (u, v) in which the linear part of T is diagonal and 
observe that the tangent vector to a n — > T(ao, cii, . . . , a n ) contains the monomials 
u n v n+1 and u n+1 v n with non-zero coefficients. However, these monomials cannot 
occur in tangent vectors to T211+1 • T(ao, . . . , a n ). This statement only involves the 
conjugation of the linear part of T by T2 n +i and is valid (with the same proof) in 
both elliptic and hyperbolic cases. Therefore (i) - (ii) are valid in both cases and 
{ao, . . • , a n } — > {bo, . . . , b n } can be inverted. □ 
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